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Note: 

The number of events quoted for DA<I>NE are based on a luminosity of 5 • 
10^^ C7n~'^s~^, which is equivalent [0] to an annual rate of 9 ■ 10^ (1.1 ■ 10^) 
tagged K"^ (Kl) (1 year = 10^ s assumed). 

Whenever we quote a branching ratio for a semileptonic K'^ decay, it stands 
for the branching ratio of the corresponding K^ decay, e.g., 

BR{K'^ -^ -n-l^v) = BR{Kl -^ -n^l^v) . 

We use the data from the Particle Data Group edition 1990 |^ throughout. 
Please contact one of the authors in case that very high precision is needed 
for a particular matrix element. We would then convert the relevant quantity 
to the newest data compilation available. 

If not stated explicitly, we always use for the low-energy constants Li, . . . , Lio 
the values displayed in table 1 in ref. 0. 

More notation is provided in appendix ^ 



1 Radiative K12 decays 

We consider the Ki2-y decay 

K+[p) ^ i+ipiMpMq) [K12,] (1.1) 

where / stands for e or fi, and 7 is a real photon with g^ = 0. Processes where 
the (virtual) photon converts into a e~^e~ or /i"'"/i~ pair are considered in the next 
subsection. The K~ mode is obtained from (|Ll|) by charge conjugation. 



1.1 Matrix elements and kinematics 

The matrix element for K~^ —>■ l^vi'-j has the structure 

T = -iGpeV:,el {FkL^ - Hn.} (1.2) 

with 

\2pq 2piq J 

l^ = u{p,)^^{l-^^)v{pi) 

W^ = {p-qY = {pi+PuY. (1.3) 

Here, e^ denotes the polarization vector of the photon with g'^e^ = 0, whereas 
y4, V stand for two Lorentz invariant amplitudes which occur in the general decom- 
position of the tensors 

F- = j dxe'^^+^'^y < I TV^Jx)i:_M I K+{p)> , I = V,A . (1.4) 

The form factor A {V) is related to the matrix element of the axial (vector) current 
in ( |1.4| ). In appendix we display the general decomposition of A'^'^, V^^ for g^ 7^ 
and provide also the link with the notation used by the PDG and in 0, 3. 

The term proportional to L^ in ( |1.2| ) does not contain unknown quantities - it 
is determined by the amplitude of the nonradiative decay K^ -^ l^vi- This part of 
the amplitude is usually referred to as "inner Bremsstrahlung (IB) contribution", 
whereas the term proportional to H^^ is called "structure dependent (SD) part" . 

The form factors are analytic functions in the complex ly^-plane cut along the 
positive real axis. The cut starts at W'^ = {Mk + 2M^Y for A (at W^ = {Mx + M^f 
for V). In our phase convention, A and V are real in the physical region of if/27 
decays, 

m] <W'^ < M|. (1.5) 



The kinematics of (spin averaged) Ki2'y decays needs two variables, for which we 
choose the conventional quantities 

X = 2pq/Ml , y = 2ppi/Ml . (1.6) 

In the K rest frame, the variable x (y) is proportional to the photon (charged lepton) 
energy, 

x = 2E^/Mk , y = 2Ei/MK , (1.7) 

and the angle ^/^ between the photon and the charged lepton is related to x and y 
by 

(1 - y/2 + A/2){1 - y/2 - A/2) _ j—^ 

1 - t//2 + A/2cos^,, '^-V^ ^"' • ^^-^^ 

In terms of these quantities, one has 

W^ = Ml(\-x) ■ (g2 = 0) . (1.9) 

We write the physical region for x and y as 

2^1 <y<\ + n 



or, equivalently, as 



where 



\-\{y^A)<x<\-]p-A) (1.10) 



< X < 1 — r; 

l-a; + -^< y <l + n (l-H) 

1 - X) 



n - mi/Mj, -< ^ g _ ^Q_2(^ ^ ^) _ (i-i2j 



1.2 Decay rates 

The partial decay rate is 



^^= 2M^ E \T\'duPsip;Pi,P.,q). (1.13) 



spins 



The Dalitz plot density 






is a Lorentz invariant function which contains V and A in the following form 



p{x,y) = PiB{x,y)+ psn{x,y)+ p,^r^{x,y) 
pm{x,y) = A,j,fjs{x,y) 



PsD{x,y) 



AsuMl 



{V + A)VsD+ (x, y) + {V- A)VsD- {x, y) 



p,^Tix,y) = A,^^MK[iV + A)f,^^+{x,y) + iV-A)f,^^-{x,y)] (1.15) 



where 



fm{x,y) 



l-y + ri 



x2 + 2(l -x)(l-n) 



x'^{x + y — 1 — ri) 
fsi>+ix,y) = [x + y-l-ri][{x + y-l){l-x)-ri] 
fsiy-{x,y) = [l-y + ri][{l-x){l-y)+ri] 
l-y + ri 



2xri{l -ri) 
X + y — 1 — ri 



fmT+{x,y) 
JiNT~ [x, y ) 
and 



x{x + y — 1 — ri) 

l-y + ri 
x{x + y — 1 — ri) 



[(1 -x){l-x -y) +ri] 
x"^ — {1 — x){l — X — y) — ri 



;i.i6) 



AjB — 

Aso = 

^INT — 

For later convenience, we note that 



4r, 



Fk_ 

Mk 



A. 



r"^ IT/ \'^n 
327r2 ^ 

4n (^ \ A. 



\M. 



A, 



a 



Sirnii-nyyFK 



K 



Mk^^ 



1.17) 



T{K -^ Ivi) . 



1.18) 



The indices IB, SD and INT stand respectively for the contribution from inner 
Bremsstrahlung, from the structure dependent part and from the interference term 
between the IB and the SD part in the amplitude. 

To get a feeling for the magnitude of the various contributions IB,SD and INT 
to the decay rate, we consider the integrated rates 



Ti= f dxdypj{x,y) ; / = SD±, INT±, IB 
Jri 



:i.l9) 



where psD = Psd+ + Psd- etc. For the region Rj we take the full phase space for 
/ ^ IB, and 

i?iB = 214.5MeV/c < pi < 231.5MeV/c . (1.20) 



Table 1.1: The quantities Xj,Nj. 80=*= and INT='= are evaluated with full phase 
space, IB with restricted kinematics ( |1.20|) . 





SD+ 


SD- 


INT+ 


INT- 


IB 




Xi 


1.67-10-2 


1.67-10-2 


-8.22 - 10-s 


3.67 -10-*^ 


3.58-10-6 


Ke2'y 


Xi 


1.18- 10-2 


1.18- 10-2 


-1.78-10-3 


1.23-10-2 


3.68-10-2 


Kf,2^ 


Ni 


2 


2 


1 


1 








for the Brenisstrahlung contribution. Here pi stands for the modulus of the lepton 
three momentum in the kaon rest system 0. We consider constant form factors V, 
A and write for the rates and for the corresponding branching ratios 



BRj 



Ni 



with 



Asu{MK{V±A)f'Xj 
Tj/T,^, = N{MKiV±A)y''Xj 



N = Asu/T,^, = 8.348 - 10" 



The values for Nj and Xj are listed in table 1.1 



:i.2ii 



:i.22) 



To estimate Tj and BR/, we note that the form factors V, A are of order 

Mk{V + A) ^ -10'^ , Mk{V - A) ^ -A ■ 10'^ . (1.23) 

From this and from the entries in the table one concludes that for the above regions 
Ri, the interference terms INT are negligible in Ke2'y, whereas they are important 
in Kfj_2-y Furthermore, IB is negligible for -^"627, because it is helicity suppressed as 



can be seen from the factor mf in A^ 



This term dominates however in K 



f,2-y 



1.3 Determination of A{W^) and ViW^) 
The decay rate contains two real functions 



^1.24) 



as the only unknowns. In Figs. ( |1.1| , |1.2| ) we display contour plots for the density 
distributions /m, . . . , /int± for I = fj,,e. These five terms have obviously very different 
Dalitz plots. Therefore, in principle, one can determine the strength of each term 
by choosing a suitable kinematical region of observation. To pin down F^, it would 
be sufficient to measure at each photon energy the interference term INT^. This 
has not yet been achieved so far, either because the contribution of INT^ is too 

^ This cut has been used in Q for it'^27, because this kinematical region is free from i^^a 
background. We apply it here for illustration also to the electron mode Ke2j- 
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Figure 1.1: Contour plots for /jb,. . . , /int± [-^^27]- The numbering on the hues 
points towards increasing modulus. The normalization is arbitrary. 



Figure 1.2: Contour plots for /jb,. . . , /int± [-^627]- The numbering on the hues 
points towards increasing modulus. The normalization is arbitrary. 
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small (in Ke2'y) , or because too few events have been collected (in K^2'y)- On the 
other hand, from a measurement of SD^ alone one can determine A, V only up to a 
fourfold ambiguity: 

SD± ^ {{V, A); -{V, A); {A, V); -{A, V)} . (1.25) 

In terms of the ratio 

IK = A/V (1.26) 

this ambiguity amounts to 

SD±^{7i^;l/7x}. (1.27) 

Therefore, in order to pin down the amplitudes A and V uniquely, one must measure 
the interference terms INT^ as well. 

1.4 Previous experiments 

K^ — > e"'"z/e7 



The PDG uses data from two experiments [^ |^ , both of which have been sensitive 
mainly to the SD+ term in ([TISl) . In [0], 56 events with E^ > 100 MeV, Ee+ > 236 

^7 



MeV and 6^+^ > 120° have been identified, whereas the later experiment 0] has 



collected 51 events with E^ > 48 MeV, E^+ > 235 MeV and 9^+^ > 140°. In these 
kinematical regions, background from K~^ -^ e'^Ueii^ is absent because i?"'"'(i^e3) = 
228 MeV. The combined result of both experiments is [| ||^ 

r(SD+)/r(K^2) = (2.4 ± 0.36) ■ 10-^ (1.28) 

For SD^, the bound 

r(SD-)/r,„,, < 1.6 ■ 10-^ (1.29) 

has been obtained from a sample of electrons with energies 220 MeV < E^ < 230 
MeV g. Using (|1.21| , |1.22|) , the result ( |T:28D leads to 



Mk\V + A\= 0.105 ± 0.008 . (1.30) 

The bound ( |1.29| ) on the other hand implies |Q 

\V -A\ / \V + A\<Vll, (1.31) 

from where one concludes |@] that 'Jk is outside the range —1.86 to —0.54, 

7k ^[-1.86,-0.54] . (1.32) 

^In all four experiments [0, §, |^, I3] discussed here and below , the form factors A and V have 
been treated as constants. 
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Table 1.2: Measured branching ratios r{K -^ l^n) /'^tot-^- The -ft'e27 data are from 



0, ^, the K^2-t data from P, ||. The last column corresponds [^ to the cut ( [1.20|) 





SD+ 


SD- 


INT+ 


SD- + INT- 


total 


Ke2-i 


(1.52 ±0.23) •10-=' 


< 1.6 -10-^ 








K^2i 


< 3 • 10-=' 




< 2.7 • IQ-'^ 
(modulus) 


< 2.6 -10-* 
(modulus) 


(3.02 ±0.10) -10-=^ 



As we already mentioned, the interference terms INT''^ m K ^ eue'y are small and 
can hardly ever be measured. As a result of this, the amplitudes A, V and the ratio 
jk determined from -ft'e27 are subject to the ambiguities ( |1.25D , ( |1.27D . 



K+ 



^J'^^^1 



Here, the interference terms INT^ are nonnegligible in appropriate regions of 
phase space (see Figs. ( [1.1| , |1.2|) ). Therefore, this decay allows one in principle to 
pin down V and A. The PDG uses data from two experiments [§, |^. In 0, the 
momentum spectrum of the muon was measured in the region (|1 . 20|) . In total 2±3.44 
SD+ events have been found with 216 MeV/c < p^, < 230 MeV/c and E^ > 100 
MeV, which leads to 

Mk\V + A |<0.16 . (1.33) 

In order to identify the effect of the SD~ terms, the region 120 MeV/c < p^ <150 
MeV/c was searched. Here, the background from i^^a decays was very serious. The 
authors found 142 K^^^ candidates and conclude that 



-1.77< Mk{V-A) <0.21. 



^1.34) 



The result ( |1.33| ) is consistent with ( |1.30| ), and the bound (|1.34|) is worse than 
the result ( |1.31| ) obtained from Ke2j- The branching ratios which follow [Q from 
( |1.33| , |1.34D are displayed in table |L2|, where we also show the ^627 results [0, Q. 
The entry SD~+INT~ for ^^27 is based on additional constraints from Ke2'y [H- 



1.5 Theory 

The amplitudes y4(Vr^) and V^(Vr^) have been worked out in the framework of various 
approaches, viz., current algebra, PCAC, resonance exchange, dispersion relations, 
.... For a rather detailed review together with an extensive list of references up to 
1976 see 0. Here, we concentrate on the predictions of V,A in the framework of 
CHPT. 



12 



1.5.1 Chiral expansion to one loop 



The amplitudes A and V have been evaluated |10|, |Tl| in the framework of CHPT 
to one loop. At leading order in the low-energy expansion, one has 



A = V = 0. 



:i.35) 



As a consequence of this, the rate is entirely given by the IB contribution at leading 
order. At the one-loop level, one finds 



A 
V 



-piLl + Ll,) 
1 1 






;i.36) 



where Lg and L^q are the renormalized low-energy couplings evaluated at the scale 
/i (the combination Lg + L\q is scale independent). The vector form factor stems 



from the Wess-Zumino term |T2| which enters the low-energy expansion at order p 
seeRef. 0. 

Remarks: 



(i) At this order in the low-energy expansion, the form factors A, V do not exhibit 
any ly^-dependence. A nontrivial l^^-dependence only occurs at the next 
order in the energy expansion (two- loop effect, see the discussion below). Note 
that the available analyses of experimental data oi K ^ lui'y decays |0) @) il Hi 
use constant form factors throughout. 

(ii) Once the combination Lg + Lio has been pinned down from other processes, 
Eq. (|1.36|) allows one to evaluate A, V unambiguously at this order in the 
low-energy expansion. Using Lg + Liq = 1.4 ■ 10~^ and F = F.,^, one has 



Mk{A + V) = 


= -0.097 


MKiV-A) = 


= -0.037 


IK = 


= 0.45 . 



;i.37) 



The result for the combination {A + V) agrees with (|1.30|) within the errors, 
while 'jK is consistent with ( p..32| ). 



We display in table [L^ the branching ratios BR/ (|1.21| ) which follow from the 
prediction (|1.37| ). These predictions satisfy of course the inequalities found from 
experimental data (see table [T^ ). 
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Table 1.3: Chiral prediction at order p^ for the branching ratios r{K -^ lvi'-))/T^, 
The cut used in the last column is given in Eq. (|1 . 20| ) . 





SD+ 


SD- 


INT+ 


INT- 


total 


Ke2i 


1.30-10-^ 


1.95-10-6 


6.64-10-1° 


-1.15 -10-s 


2.34 - 10-6 


K^2j 


9.24 ■ 10-^ 


1.38-10-6 


1.44-10-5 


-3.83-10-5 


3.08 - 10-3 



1.5.2 ly ^-dependence of the form factors 

The chiral prediction gives constant form factors at order p"^. Terms of order p^ 
have not yet been calculated. They would, however, generate a nontrivial W'^ - de- 
pendence both in V and A. In order to estimate the magnitude of these corrections, 
we consider one class of p^ - contributions: terms which are generated by vector and 
axial- vector resonance exchange with strangeness 
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V{W^' 



V 



1 - W^Mk^ 



A{W^ 



wyMK,' 



:i.38) 



where V,A are given in ( |1.36D . We now examine the effect of the denominators in 
( |1.38|) in the region y > 0.95, x > 0.2 which has been explored in K~^ -^ e^Vel 0- 
We put rrie = and evaluate the rate 



dPjx) 
dx 



N.. 



y=0.95 



Pso+{x,y)dy 



:i.39) 



where iVt^t denotes the total number of K~^ decays considered, and T^J = 1.24- 10 



sec. 



dP(x) 



The function J-^' is displayed in Fig. ( |1.3[ ) for three different values of M^* 
and Mkj^, with N^^^ = 9 - 10^. The total number of events 



Np= f dP{x) 

Jx=0.2 



[lAO) 



is also indicated in each case. The difference between the dashed and the dotted line 
shows that the nearby singularity in the anomaly form factor influences the decay 
rate substantially at low photon energies. The effect disappears at a; — >■ 1, where 



W 



Ml{l 



X) 



0. To minimize the effect of resonance exchange, the large-x 



region should thus be considered. The low-x region, on the other hand, may be used 
to explore the py^-dependence of V and of A. For a rather exhaustive discussion of 
the relevance of this W^ - dependence for the analysis of Ki2^ decays we refer the 
reader to Ref. HI. 
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Figure 1.3: The rate dP{x)/dx in (|1.39|) , evaluated with the form factors ( |1.38|) and 
N,,, = 9 • 10^ The sohd hue corresponds to Mk* = 890 MeV, Mk, = 1.3 GeV. 
The dashed hne is evaluated with M^* = 890 MeV, Mk^ = C)0 and the dotted line 
corresponds to Mx* = M^i = C)0. The total number of events is also indicated in 
each case. 



1.6 Comment on tensor couplings 



Bolotov et al. WM have analyzed radiative pion decays vr ^ e Uf.'y in flight (~ 80 
events) in a wider kinematical region than was explored in the high-statistics exper- 



iment of Bay et al. [^ (where ~ 700 



IT 



e'^Vel events had been observed). The 



theoretical branching ratio, calculated with the standard V — A coupling, differs 
from the measured one by more than three standard deviations. This discrepancy 
may be avoided by adding to the standard matrix element the amplitude of a ten- 
sorial interaction [|l^. Belyaev and Kogan |jT^ and Voloshin [|l8l have pointed out, 
however, that in the standard model the induced tensor coupling is too small to 
generate the rate observed in ref. |14 . 



Gabrielli 19 



has worked out the effect of tensor couplings for K^ -^ I i^n 
decays. Using the above quoted values for the form factors A and V and a tensor 
coupling of a size suggested to explain the data in Ref. |T^, he finds a < 30% effect 
in the partial decay rates (the exact size depends on the chosen coupling, channel, 
decay region,...). The author then suggests that these effects may be accessible to 
detection at high precision experiments carried out at DA$NE. 

We wish to point out that this may be difficult for the following reason. The 
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calculation of the decays K^ -^ l^^n presented in this section is based on the 
one-loop formulae for the decay matrix elements. Higher-order effects may well be 
sizeable, see e.g. figure [L^. There, it is explicitly seen that the effect of resonance 
exchange is ~ 30% in particular regions of phase space. Therefore, in order to 
identify effects due to tensor couplings, one first has to pin down the contribution 
from higher-order effects in CHPT. This is not an easy task to achieve to the accuracy 
required. On the other hand, it is of course needless to say that the finding of a 



tensorial coupling of the size suggested in Ref. |]T6i would be spectacular 



1.7 Improvements at DA$NE 

Previous experiments have used various cuts in phase space in order (i) to identify 
the individual contributions IB, SD^*", INT^ as far as possible, and (ii) to reduce the 
background from Ki^ decays. This background has in fact forced so severe cuts that 
only the upper end of the lepton spectrum remained. 

The experimental possibilities to reduce background from Ki^^ decays are pre- 
sumably more favourable with today's techniques. Furthermore, the annual yield of 
9 ■ 10^/C"'" decays at DA$NE is more than two orders of magnitude higher than the 
samples which were available in 0, ^ |^, |^ . This allows for a big improvement in the 
determination of the amplitudes A and V , in particular in K^2'y decays. It would be 
very interesting to pin down the combination Lg + Lio of the low-energy constants 
which occur in the chiral representation of the amplitude A and to investigate the 
VT^-dependence of the form factors. 
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2 The decays K^ -^ l^ul'+l'- 

Here we consider decays where the photon turns into a lepton-anti-lepton pair, 



-+ 


-> e~^z/yU^/i" 


-+ 


-> fi'^ue^e' 


-+ 


-^ e^ve^e~ 


-+ 


-^ /i^Z/;U+/i 



(2.1) 
(2.2) 
(2.3) 
(2.4) 



2.1 Matrix elements 

We start with the processes (|2.1| ) and 

(/,/') = (e,/i)or(/i,e). (2.5) 

The matrix element is 

T = -iGpeV:,!, {FkV - H'^'l^} (2.6) 

where 



2pg — g2 2pig + g^ 



miu(p^)(l + 75) <^ - —^ } v{pi 



with 



-A2{q^gP^ - qfq^) - A^^qWqP - q^WP)^ (2.7) 

The form factors Aj(g^, VT^), V^i(5'^, W^^) are the ones defined in appendix |y. Fy{q^) 
is the electromagnetic form factor of the K^ . Finally the quantity e^ stands for 

e^ = ^u{p2)l^v{pr) , (2.9) 

and the four-momenta are 

q = P1+P2, W = pi+pi,=p-q (2.10) 

such that q^e'^ = 0. 

In order to obtain the matrix element for ( ^.3p and ( |2.4| ), 

i^+(p) ^ /+(pO^(p.)/''(Pi)r(p2) , (2.11) 
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one identifies mi and m'l in ( p.6| ) and subtracts the contribution obtained from 
interclianging pi ^^ pi : 

(PuPi) -^ {PhPl) 
q -> P1+P2 
W -> p-q = p^+p^ . (2.12) 



2.2 Decay distributions 

Tlie decay widtli is given by 

'^^ = o/i#_.^o^A8 Y. \T^\dLiPs{p;PhPu,Pi,P2) (2.13) 



2MK{27Ty 



spms 



and tlie total rate is tlie integral over tliis for the case / 7^ I'. For the case I = I' 
the integral has to be divided by the factor 2 for two identical particles in the final 

state. 

We first consider the case where / 7^ /' and introduce the dimensionless variables 

_ 2pq_ _'2pip _ _g^ _^ /_^ /^..x 

"'"Ml-' ^" Mi' '^" Mi' '■'"Mi''"' "Mi • ^^-'^^ 



Then one obtains, after integrating over pi and p2 at fixed q [EO 



dV K+^i+vi'+v- = a'^Gp\Vus\'^Ml^F{z,r'i) < - ^ T*T'' > dxdydz 



spins 



r = M^2|^^;^M_;^M^^^|_ ^2.15) 

This result allows one to evaluate, e.g., the distribution dT/dz of produced I'^l'^ 
pairs rather easily. The kinematically allowed region is 



(2.16) 
with 





4r[< 


z <l+ri- 


2v^ 




2v^< 


X < 1+z- 


ri 




A-B< 


y < A + B 




A 


(2- 


x){l + z + ri- 


-X) 


(1 + 


2{l + z-x) 
z-x-ri)\fx 




p 


2-4z 






2(l + z-x) 
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(2.17) 



The case / = /' is slightly more elaborate. We feel that it does not make sense to 
display the term Jospins I^P because it is of considerable complexity in the general 
case when all the form factors Ai, V\ and Fy are g^ and W'^ dependent. The 
expression together with the Monte Carlo program to do the phase space integrals 
is available on request from the authors. 

2.3 Theory 

The form factors Aj, V\ and Fy have been discussed in all kinds of models, Vector 
Meson Dominance, hard meson, etc.. For a discussion see Ref. 0|. We will restrict 
ourselves to the predictions in the framework of CHPT. 

To leading order we have 

Fi = 

A^ = A2 = A3 =0 . (2.18) 

We also have Fy = 1. The rate here is entirely given by the inner Bremsstrahlung 



contribution. At the one-loop level several form factors get non-zero values |TI 

1 



Vi 



Stt^F 



Ai = -j{Ll + L\,) 



Ao 



2FK{F^{q')-l) 
q2 



A3 = 
F^(g2) = l + HUq^) + 2HKK{q^) . (2.19) 

These results obey the current algebra relation of Ref. 0. The function Fy[q^) 
does, however, deviate somewhat from the linear parametrization often used. The 
function H{t) is defined in appendix |B[ 

The fact that the form factors at next-to- leading order could be written in terms 
of the kaon electromagnetic form factor in a simple way is not true anymore at the 
p^ level. The Lagrangian at order p^ contains a term of the form 

ti [d^FI^U^D^Frp.u] (2.20) 

that contributes to A2 and ^3 but not to the kaon electromagnetic form factor, 

2.4 Numerical results 

We have calculated the rates for a few cuts, including those given in the literature. 
For the case of unequal leptons, the results are given in table |2]l| for the decay K^ —^ 
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Table 2.1: Theoretical values for the branching ratios for the decay K^ 
for various cuts. 



fi^ue^e 





tree level 


form factors as given by CHPT 


full phase space 


2.49 ■ 10-5 


2.49 - 10-5 


z < 10"^ 


2.07-10-5 


2.07-10-5 


z > 10-=^ 


4.12-10-6 


4.20 - 10-6 


z > (20 McV/MkY 


3.15 -lO-'^ 


3.23 ■ 10-6 


2 > (140 McV/MkY 


4.98 -10-s 


8.51 - 10-^ 


X > 40 McV/Mk 


1.58-10-5 


1.58-10-5 



/i+z/e+e-. These include the cuts used in Refs. |^ and [^, x > 40 McV/Mk and 
z > (140 MeV/Mx)"^, respectively. It can be seen that for this decay most of the 
branching ratio is generated at very low electron-positron invariant masses. As can 
be seen from the result for the cuts used in Ref. lETI, the effect of the structure 



dependent terms is most visible at high invariant electron-positron invariant mass. 
Our calculation, including the effect of the form factors agrees well with their data. 
We disagree, however, with the numerical result obtained by Ref. pO| by about an 
order of magnitude. 

For the decay K~^ — > e^ufi^fi" , we obtain for the tree level or IB contribution a 
branching ratio 

BRib{K^ -^ e+z/^+/i-) = 3.06 - 10-^^ (2.21) 



and, including the form factors, 

BRtotal (K 



e vfi fi 



1.12- 10" 



(2.22) 



Here the structure dependent terms are the leading-contribution since the inner 
Bremsstrahlung contribution is helicity suppressed as can be seen from the factor 
mi in L^. 

For the decays with identical leptons we obtain for the muon case a branching 
ratio of 

BRtotai{K+ -^ /i+z//i+;u-) = 1.35 - 10-^ (2.23) 

for the full phase space including the effects of the form factors. The inner Brems- 
strahlung or the tree level branching ratio for this decay is 



BRib{K+ -^ fi+ufi+^-) = 3.79 - 10-^ 



(2.24) 



For the decay with two positrons and one electron the integration over full phase 
space for the tree level results is very sensitive to the behaviour for small pair masses. 
We have given the tree level and the full prediction, including form factor effects in 
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Table 2.2: Theoretical values for the branching ratios for the decay K^ -^ e+z/e+e 
for various cuts. 





tree level 


form factors as given by CHPT 


full phase space 


?«4-10-9 


1.8-10"^ 


z, zi > 10"^ 


3.0-10-^" 


1.22 -10"'^ 


z, zi > (50 MeV/MK^ 


5.2-10-11 


8.88 - 10-s 


z, zi > (140 MeV/MK^ 


2.1-10-12 


3.39 - 10-s 



table 2.2. The cuts are always on both invariant masses 



z 



{p,+P2f/Ml 



{pi + p^y/M^ 



K 



The values for the masses used are those of K^ and vr" 
the values given in table 1 in Ref. P|, 



(2.25) 
For Lg and Lio we used 






6.9-10 



-3 



-5.5-10 



-3 



(2.26) 



2.5 Present experimental status 

Only decays with an electron positron pair in the final state, decays ( p.2| ) and 
have been observed. 

Both have been measured in the same experiment [^. The decay K^ — >■ 
/i+z/e+e- was measured with a branching ratio of (1.23 ± 0.32) - 10"'^ with a lower 
cut on the electron positron invariant mass of 140 MeV. The measurement is com- 
patible with our calculation including the form factor effects for the relevant region 
of phase space. This measurement was then extrapolated |^ using the result of ||2^ 
to the full phase space. Since we disagree with that calculation, we also disagree 
with the extrapolation. 

In the same experiment, 4 events of the type K'^ -^ e^ve^e^ were observed where 
both electron positron pair invariant masses were above 140 MeV . This corresponds 
to a branching ratio for this region of phase space of (2.8l;^;4) ' 10-*^. This result 
is compatible within errors with our calculation, see table |2.2| . The matrix element 
of Ref. [^] was again used for the extrapolation to full phase space[^]. Apart 
from our numerical disagreement, the calculation of Ref. pO|] was for the case of 
non-identical leptons and cannot be applied here. 

For the decay K^ -^ /i+z//x+/x~ an upper limit of 4.1 - 10~^ exists [B^. This upper 



limit is compatible with our theoretical result, Eq. (|2.23|) . 
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The decay K^ -^ e^u^^^^ has not been looked for so far and should be within 
the capabilities of DA$NE given the branching ratio predicted in the previous sub- 
section. This decay proceeds almost entirely through the structure dependent terms 
and is as such a good test of our calculation. 

2.6 Improvements at DA$NE 

The decays discussed in this subsection, K^ -^ l^vl^l~ , are complementary to the 
decays K^ -^ /+Z/7. As was the case for the analogous decay, vr"'' -^ e^ue^e^ |l23| , 
it may be possible to explore phase space more easily with this process than with 
K^ — > l~^i"j to resolve ambiguities in the form factors. 



As can be seen from our predictions, tables |2.1| and p.2| , all the decays considered 
in this subsection should be observable at DA$NE. Large improvements in statistics 
are possible since less severe cuts than those used in the past experiments should be 
possible. In the decays with a yU^/i^ pair and the decay K~^ —^ e^ve^e~ the effects 
of the form factors are already large in the total rates and should be easily visible at 
DA<I>NE. In the decay K^ — ^ ix^ue^e~ most of the total rate is for small invariant 
mass of the pair and is given by the inner Bremsstrahlung contribution. There are, 
however, regions of phase space where the form factor effects are large and DA$NE 
should have enough statistics to be able to study these regions. 
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3 if/3 decays 

The decay channels considered in this subsection are 

K^p) -. Ap'riPiMP.) [K,",] (3.1) 

K\p) -> n-ip'ripiMp.) [K^s] (3.2) 

and their charge conjugate modes. The symbol / stands for /i or e. We do not 
consider electromagnetic corrections and correspondingly set a = throughout this 
subsection. 

3.1 Matrix elements and kinematics 

The matrix element for K^ has the general structure 

T = ^V:rF,^p',p) (3.3) 



with 



l^ = u{p,)i^{l-j,)v{pi) 
F,Hp,p) = <n'ip')\V;~^'iO)\K^ip)> 

= ^[(p'+p)./r^"(t)+(p-p')./5"^"(t)]. (3.4) 

To obtain the matrix element for K^^ , one replaces F^ by 

F,'ip\p) = <n-ip')\V;-^\0)\K'ip)> 

= {p+p),ff^~it) + {p-p'),f^'--{t). (3.5) 

The processes (|3.1|) and ( |3.2| ) thus involve the four Ki^ form factors f^ '" (t), 
f^ "^ [t) which depend on 

t = {p-py = {pi+Pu)', (3.6) 

the square of the four momentum transfer to the leptons. 

Let f^'^ = f^ '^ or f^ '^ . f^^ is referred to as the vector form factor, because 
it specifies the P-wave projection of the crossed channel matrix elements < | 
V^~*^(0) I if+,7r° in >. The S-wave projection is described by the scalar form 
factor 

ft it) = ft it) + ^^^^, /_^-(t) . (3.7) 



Analyses of Ki^ data frequently assume a linear dependence 

t 



/f,o(i) = frio) 
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1 + A, 



,0 






(3. 



For a discussion of the validity of this approximation see [^ |^ and references cited 
therein. Eq. ( p.8|) leads to a constant f'^'^if) , 

/^'^(i) = /^^(O) = /r(0)(Ao - A+) ^, - ■ (3.9) 

The form factors f±Q{t) are analytic functions in the complex t-plane cut along the 
positive real axis. The cut starts at t = {Mk + M^)^. In our phase convention, the 
form factors are real in the physical region 

m] <t< (Mk - M^f. (3.10) 

The kinematics of (spin averaged) Ki^ decays needs two variables, for which we 
choose 

y = 2ppi/Ml , z = 2pp'/Ml = {-t + Ml + MD/M], . (3.11) 

In the K rest frame, y {z) is proportional to the charged lepton (pion) energy, 

y = 2Ei/MK , z = 2EJMk • (3.12) 

The physical region for y and z is 

2^/rl< y < 1 +n -r^ 
A{y)-B{y)< z < A{y) + B{y) 

A{y) = { 2-y){l + n + r^-y)/\2{l + ri-y)] 

B{y) = ^y^ - 4n(l + n ~ r^ - y)/[2{l + n - y)] 

ri = mi /Ml,T^ = Ml/Ml. (3.13) 

or, equivalently, 

2^ < z <l + r^-ri 
C{z)-D{z)< y <C{z) + D{z) 

C{z) = {2-z){l + r^ + ri-z)/[2{l+r^-z)] 

D{z) = ^z^ - 4r^(l +r^-ri- z)/[2{l + r^ - z)] . 

(3.14) 

3.2 Decay rates 

The differential decay rate for K^ is given by 

^^ = o?i# to \5 ^ \T\'^dLiPs{p;PhPu,p') ■ (3.15) 

2Mx(27r)& spin. 
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The Dalitz plot density 

0(11 7.^ = 

dydz 2567r3 



pM = ^ = ^,^\T\' (3.16) 



spins 



is a Lorentz invariant function which contains f^ ^ in the following form, 



M%Gl I Vus 



|2 



pii'+7r0\2 , tdcK+tt''cK+-k° , rif-eK+T:°\2 



P(y^ ^) = "2567r3 [^(/+ ) + Bf^--"f^--" + Cif^-^'^l (3.17) 



with 



A{y,z) = A{z + y-l){l-y)+ri[Ay + 3z-3]-Ar^ + ri{r^-ri) 

B{y,z) = 2ri{3 - 2y - z + ri - r^) 

Ciy,z) = ri{l + r^-z-ri). (3.18) 

The quantities {A, B, C) are related to the ones quoted by the PDG [^ by 

{A,B,C) = ^{A,B,CUo . (3.19) 



To obtain the rate for i^^g, one replaces in ( |3.17| ) f^ ^ by v^/±' "^ 



For convenience we also display the K^^^/ Kf,^ rates evaluated in the approxima- 
tion (B.SI) for the form factors, 



+ + ^ 0.645 + 2.087A+ + 1.464Ao + 3.375A^ + 2.573Ag 

U^^sJ/U^eaJ - 1 + 3.457A+ + 4.783A2+ 

. 0-645 + 2.086A, + 1.459Ao + 3.369AU2.560Ag . (3.20) 
^ ^'^" ^ "^^ 1 + 3.456A+ + 4.776A^ ^ ' 

We have used the physical masses [^ in evaluating these ratios and Mt,+ to scale 
the slope in both cases. The terms linear and quadratic in Aq are proportional to mf 
and therefore strongly suppressed in the electron case. We do not include them in 
the denominators, because these coefficients are smaller than 10"^. The interference 
term AoA+ is absent by angular momentum conservation. Furthermore, one has 

d,d.Afe..) = {«;«^^3 Igj (3,21) 

3.3 Determination of the K12, form factors 



Measurements of the Dalitz plot distribution ( p.l7|) of K^^^ data allow one in principle 
to pin down the form factors (up to a sign) in the range im? < t < {Mk — M^^^"^ 



Measuring the K^^/Kes branching ratio and then using (|3.2CI| ) gives a relationship 
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between A+ and Aq which is vahd in the approximation ( |3.8| ). Furthermore, muon 
polarization experiments measure the weighted average of the ratio f^^{t)/ f^'^{t) 
over the t range of the experiment [0, p5[ . On the other hand, the electron modes 



i^es are sensitive to f^"^ only, because the other contributions are suppressed by the 
factor {mjMKf ^ 10-^ see eqs. ( pTzD , ( |3:T8| ). 

Isospin breaking effects in f^ '" (0) and f^ '^ (0) play a central role in the 
determination of the Kobayashi-Maskawa matrix element Vus from i^es data, see 



| 26[| for a detailed discussion of this point. In the following we concentrate on the 
measurement of the slopes A+^o- 



3.4 Previous measureraents 

We refer the reader to the 1982 version of the PDG |2^ for a critical discussion 
of the wealth of experimental information on A:^q. Here we content ourselves with 
a short summary. 

A'es-experiments 

The A+ values obtained are fairly consistent. The average values are 



K:, ■■ Ah 






A, 



0.028 ±0.004 Ref.[3] 
0.030 ±0.0016 Ref.iHl 



(3.22) 



ir^3-experiments 

The result by Donaldson et al. [^| 



A+ = 0.030 ±0.003 
Ao = 0.019 ±0.004 



(3.23) 



dominates the statistics in the K^.^ case. 



The A 4- value 



is consistent with 



the Kez value ( |3.22| ). However, the situation concerning the slope Aq is rather 
unsatisfactory, as the following (chronological) list illustrates. 



0.0341 

0.050 

0.039 

0.047 

0.025 

0.019 



± 
± 
± 
± 
± 
± 



0.0067 

0.008 

0.010 

0.009 

0.019 

0.004 



[29] 

m 
m 



m 
m 



(3.24) 



The x^ fit to the K^ data yields A+ = 0.034 ± 0.005, Aq = 0.025 ± 0.006 with a 



X^ I DF = 88/16 1^, p. 76]! The situation in the charged mode K^^ is slightly better 



27 



^Please note that the most recent measurements of A+_o go back to 1981 ||]! 
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3.5 Theory 

The theoretical prediction of Ki^ form factors has a long history, starting in the 
sixties with the current algebra evaluation of f^ ^ . For an early review of the 
subject and for references to work prior to CHPT evaluations of f± we refer the 
reader to [^ (see also Ref.[^). Here we concentrate on the evaluation of the form 
factors in the framework of CHPT. We restrict our consideration to the isospin 
symmetry limit rriu = ma, as a result of which one has 

/S-"(t) = f^:-\t) ^ /±,o(t) ■,m^ = m,. (3.25) 



3.5.1 Chiral prediction at one- loop order 



In Ref. [^, the vector current matrix elements < M' \ q'j^^q \ M > have been 
calculated up to and including terms of order t = {p' — pY and of order rriu, rrid 
and rUs in the invariant form factors. For reasons which will become evident below, 
we consider here, in addition to the Ki^ form factors, also the electromagnetic form 
factor of the pion 

< vr+(p') I V^^JO) I 7r+(p) >= ip' + pYF^it). (3.26) 

The low-energy representation for Fy{t) |2^, |3^ and /+(t) ||2^ reads 



U{t) = l + ^HKAt) + lHKr,it). (3.27) 

The quantity H{t) is a loop function displayed in appendix |B[ It contains the low- 
energy constant Lg. The indices attached to H{t) denote the masses running in the 
loop. 

Since Lg is the only unknown occurring in Fy{t) and in f+{t), we need ex- 
perimental information on the slope of one of these two form factors to obtain a 
parameter-free low-energy representation of the other. 

The analogous low-energy representation of the scalar form factor is 

1 / A^ \ - 

/oW = l + ^(5t-2EK,-3^) J^,(t) 

+ ^ , ^r> I 3t — 222xtt ' 

(3.28) 



1 + 


1 / 

8F2 1^ 


1 


(o. 


24F2 \^' 


t 


(Fk 


^K. V F^ 
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Figure 3.1: The vector and scalar form factors /+(t) and /o(t). 
The function J(t) is hsted in appendix |B|, and T^xn and A^tt stand for 

AKn = Ml -Ml . (3.29) 



The measured value |26] Fx/ F.,^ = 1.22±0.01 may be used to obtain a parameter-free 



prediction of the scalar form factor /o(t). 

3.5.2 Momentum dependence of the vector form factor 

In the spacelike interval \f—i < 350 MeV the low-energy representation ( p.27| ) 



for the electromagnetic form factor Fy(t) is very well approximated by the first two 
terms in the Taylor series expansion around t = 0, 

F^(t) = l + i<r2>- t + ... . (3.30) 



Likewise, the linear approximation 

/+W = /+(0){n-i<r2>^'^t + --j (3.31) 



reproduces the low-energy representation (|3.27|) very well, see Fig. p.l| . This is in 
agreement with the observed Dalitz plot distribution, which is consistent with a 
form factor linear in t. The charge radii are 

9 12^Q 1 { , Ml ^ Ml 1 
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< r^ >--< r- >- ^— Uh, i^\ + 3/.. ("^^ 






+ n^^T^ + 7,i^Tjt-^\ (3-32) 

where 



1 (x3 - 3x^ - 3x + 1) , l/x + l\2 1 

^^(") = 2 ^^Tp ^""+2l^ "3 • ^'-''^ 

To evaluate these relations numerically, we use the measured charge radius of the 
pion: 

< r^ >l= 0.439 ± O.OOSW [37] (3.34) 

as input and obtain the prediction 

A+ = -M^+ < r^ >^^= 0.031 (3.35) 

6 



in agreement with the experimental results (|3.22|) , ( p.23| ) . From this (and from 



the considerably more detailed discussion in Ref. l2J]), one concludes, in agreement 
with other theoretical investigations [^, that the measured charge radii < r^ >^ 
and < r^ >y'" are consistent with the low-energy prediction. 



3.5.3 Momentum dependence of the scalar form factor. 
Dashen-Weinstein and Callan-Treiman relations 

In the physical region of Ki^ decay the low-energy representation ( |3.28| ) for the 
scalar form factor is approximated by the linear formula 

/o(t) = /+(0) |l + 1 < r^ >f- t + . . j (3.36) 



to within an accuracy of 1 %. (See Fig. |3.1|) . The curvature generated by higher- 
order terms is also expected to be negligible in the physical region of the decay 
For the slope < r^ >§'" one obtains 



<^'>r - MF^(t-i)-^*=^o(*-"-) 



Ml] 19Mf, + 3M^ ^ M^ 



(3.37) 



'^ We do not quote an error for the result ( |3.35 ), because one should estimate higher-order chiral 
corrections for this purpose. 
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Figure 3.2: The normalized slopes of the vector and the scalar form factors. Curve 
1: the normalized slope M^+(i/+(t) /(it. Curve 2: the normalized slope M^+(i/o(t) /(it. 
Near the ttK threshold to = (M^ + M^)^, the vector form factor behaves as /+(t) = 
/+(to) + 0[{t — to)], whereas /o(t) = /o(to) + 0[{y/t — to)]. The slope of the scalar 
form factor is therefore singular at t = {Mk + M^)^. 



where 



3 /l + xX^ 3x(l + x) 
2 U-xj ^ fl-x) 



^2(3;) = -{- ::) + ,, X3 ln^> 



/12(X) = /l2 (^) , /^2(1) = 1, 

m = (mu + mrf)/2 . (3.38) 

This (parameter-free) prediction is a modified version of the Dashen-Weinstein re- 
lation [^, which results if the nonanalytic contribution 62 is dropped. Dashen, Li, 
Pagels and Weinstein [^] were the first to point out that the low-energy singularities 
generated by the Goldstone bosons affect this relation. The modified relation is for- 
mulated as a prediction for the slope of /o(t) at the unphysical point ti = Mj^ + M^. 
Their expression for this slope however has two shortcomings: (i) it does not account 
for all corrections of order Ai; (ii) The slope at ti differs substantially from the slope 
in the physical region of the decay ||2^, ^ , see Fig. |3^ . 

Algebraically, the correction 62 is of the same order in the low-energy expansion 
as the term involving Fk/Fj^ — 1. Numerically, the correction is however small: 62 
reduces the prediction by 11 %. With Fk/Fj^ = 1.22 ± 0.01 the low-energy theorem 
( p.37| ) implies 



< r^ >f ^ = 0.20 ± O.OSfm^ 
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Ao = -M2+ <r2 >f^= 0.017 ±0.004 (3.39) 

6 

where the error is an estimate of the uncertainties due to higher-order contributions. 
The prediction ( |3.39| ) is in agreement with the high-statistics experiment [^ quoted 
in (|3.23| ) but in flat disagreement with some of the more recent data hsted in ( p.24|) . 



In the formulation of Dashen and Weinstein 39 , the Callan-Treiman relation 



1^2] states that the scalar form factor evaluated at t = M^ — M^ differs from Fk/F^^ 
only by terms of order m^, 771^: the quantity 

Act = foiMl -Ml)-^ (3.40) 



is of order m. Indeed, the low-energy representation ( p.2S ) leads to 

Act = -^ \Jk.{MI - Ml) + -Jk^{MI - Ml)^ + 0{rhms) . (3.41) 

Numerically, Act = —3.5 ■ 10^'^. The Callan-Treiman relation should therefore 
hold to a very high degree of accuracy. If the form factor is linear from t = to 
t = M]^ — Ml then the slope must be very close to 

in agreement with ( p.39| ) and with the experimental result of Ref. |Q, but in 
disagreement with, e.g., the value found in Ref. [^. We see no way to reconcile 
the value Aq = 0.050 with chiral symmetry. 



3.6 Comment on tensor couplings 



S.A. Akimenko et al. [Q have investigated the general form of the matrix element 
for Ke3 decays, obtained by adding scalar and tensor couplings to the standard V —A 
interaction. Analyzing the Dalitz-plot distribution of 3.2 ■ 10^ K^ -^ ii^e^Ve events, 
they find that the presence of scalar and tensor couplings or nonlinearities in the 
form factor /+ cannot be excluded. DA<I>NE may be an ideal place to check this 
claim (~ 4 ■ 10* K^ — ^ ir^e^Ve events in one year). However, the same proviso as 



in the case of radiative K12 decays should be made (see subsection pTq) : before any 
firm conclusion can be drawn, one has to estimate the effect of higher-order terms 
in the chiral calculation. In the present case this may be less difficult to achieve 
than in radiative kaon decays, as only one form factor comes into play, which, in 
addition, depends on one kinematical variable only. 
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Table 3.1: Rates of if;3 decays. The number of events in the third column corre- 
sponds to those data which are of relevance for the determination of the slope Aq of 
the scalar form factor. 





y events 






branching 
ratio 


Particle Data 
Group 


DA$NE 

1 year 


improve- 
ment 


K+ -^ ttV+^^m 


3.18-10-2 


10^ 


3-108 


3-103 


Kl —* TT^/i^Z/ 


27 ■ 10-2 


4-10'' 


3-10^ 


70 



3.7 Improvements at DA$NE 

DA$NE provides the opportunity to improve our knowledge of K13 decays in a 
very substantial manner - in particular, it should be possible to clarify the issue 
of the slope Aq of the scalar form factor /q. To illustrate, we compare in table 



3TT| the hitherto obtained number of events (third column) with the expected ones 
at DA$NE (fourth column). The last column displays the remarkable increase in 
statistics obtainable at DA$NE. 
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4 Radiative i^/3 decays 

The decay channels considered in this subsection are 

K\p) -. n-ip'r{p,Mp^Uq) [irg 

and the charge conjugate modes. We only consider real photons (g^ = 0). 

4.1 Matrix elements 

The matrix element for K^ has the general structure 

T = ^eV:,e'^iqr!^iV;^-Apuip^)Yil-lMPi) (4.1) 

'^PiQ J 



The diagram of Fig. ^]T|.a corresponding to the first part of Eq. ( |4.1| ) includes 
Bremsstrahlung off the K~^. The lepton Bremsstrahlung diagram of Fig. [4.1| .b is 
represented by the second part of Eq. ( [4.1| ). The hadronic tensors V^^.A^^ are 
defined as 

J+ =ijd'xe^^%7r\p') I T{V;-(x)Jr^(0)} I ir+(p)), / = V, A (4.2) 

F^ is the K^ matrix element 

F^ = {n\p')\Vt'\(^)\K+{p)). (4.3) 

The tensors V^"^ and A^^ satisfy the Ward identities 

q'v^. = K (4.4) 

q'A^, = 

leading in turn to 

q'K = ' (4-5) 

as is required by gauge invariance. 

For Kj^^ , one obtains the corresponding amplitudes and hadronic tensors by 
making the replacements 

K+ -^ K^, TfO^TT- 

V,X - V^^, A;,^A% (4.6) 

F+ ^ po A+ -^ 4° 
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To make the infrared behaviour transparent, it is convenient to separate the 

^+ T/O 



tensors V"^, V° into two parts: 



V+ = F+ + i-^F+ (4.7) 



" IJ.U " fiu ^ , -•- u- 

Due to Low's theorem, the amphtudes V^j^ are finite for g ^ 0. The axial amph- 
tudes At^^ are automatically infrared finite. The Ward identity ( [4.4|) implies that 

the vector amplitudes V^j^ are transverse: 

q^V+f = 0. (4.8) 

For on-shell photons, Lorentz and parity invariance together with gauge invari- 
ance allow the general decomposition (dropping the superscripts +,0 and terms that 
vanish upon contraction with the photon polarization vector) 

+^3 {p',W^ - ^W,w)j + V, l^p'^pl - ^W,p'}j (4.9) 

V = ^^M-pal^iP'"?'' + v42gW^) + te^^p^p'^q^WiA^W, + A^p'J 
Fu = Cip', + C2{p-p)u 

W = Pl+Pu- 

With the decomposition ( f4.7| ) we can write the matrix element for Kjf^ in ( [4.1| ) in 
a form analogous to Eq. ( |1.2|) for Ki2'y: 



T = %eK>^(g)*|(V;t-A+)ll(p.)7^(l-75)t;(p/) (4.10) 



+ Ku[Pu)i (1 - 75) 



pq 2piq 



v{pi 



The four invariant vector amplitudes Vi, . . . , V4 and the four axial amplitudes 
Ai,...,Ai are functions of three scalar variables. A convenient choice for these 
variables is 

E^ = pq/MK, E^ = pp'/Mk, W = Vw^ (4.11) 

where W is the invariant mass of the lepton pair. The amplitudes Ci,C2 can be 
expressed in terms of the Ki^ form factors and depend only on the variable {p—p'Y = 
M"^ + M^ — 2MkEt^. For the full kinematics of Ki^^ two more variables are needed, 
e.g. 

El = ppi/Mk, X = piq/Ml. (4.12) 
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Figure 4.1: Diagrammatic representation of the Kt amplitude. 

The variable x is related to the angle di^ between the photon and the charged lepton 
in the K rest frame: 



xMl = E^{Ei - ^Ef - mf cosOi^). (4.13) 

T invariance implies that the vector amplitudes \4,...,V4, the axial amplitudes 
Ai, . . . ,A4 and the Ki^ form factors 6*1,6*2 are (separately) relatively real in the 
physical region. We choose the standard phase convention in which all amplitudes 
are real. 

For 61^ — > (coUinear lepton and photon), there is a lepton mass singularity in 
( [4.1|) which is numerically relevant for I = e. The region of small E^, 61^ is dominated 
by the Ki^ matrix elements. The new theoretical information of Ki^^ decays resides 
in the tensor amplitudes V^y and A^^. The relative importance of these contributions 
can be enhanced by cutting away the region of low E^, 61^. It may turn out to be 
of advantage to reduce the statistics by applying more severe cuts than necessary 
from a purely experimental point of view. 



4.2 Decay rates 

The total decay rate is given by 

T{K^nlu^) = —-l—^JduPs{p;p',Pi,Pu,q)T.\T\' (4.14) 



spins 



in terms of the amplitude T in ( |4.1| ). The square of the matrix element, summed 
over photon and lepton polarizations, is a bilinear form in the invariant amplitudes 
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Vi, . . . , V4, Ai, . . . , A4, Ci, C2. Pulling out common factors, we write ( ^.14| ) in the 
form 

T{K -. nlu^) = ,J,|,7' / duPs{p;p',Pi,P.,q) SM , (4.15) 



{2ixyMK J 

where SM is the reduced matrix element. For the actual numerical calculations, we 
have found it useful to employ a tensor decomposition different from the one in Eqs. 
Q and §3) 



V, 



flU 



+B,W,W, + Bep'^W, + Brp'^pl . 



(4.16) 



One advantage is that ( |4.16| ) applies equally well to both charge modes while ( [4.7|) 
does not. In the numerical evaluation of the amplitudes, gauge invariance can of 
course be used to express three of the Bi in terms of the remaining ones and of 
Ci, C'2. 

To get some feeling for the magnitude of the various decay rates, let us first con- 
sider the tree level amplitudes to lowest order p^ in CHPT. With the sign conventions 



of Ref. [0, these amplitudes are |Tl|, ^| 



K, 



i37 • 



V. 



Hv 



A 



flU 



1 

71 



1 

V2 



9^,u + 



{p' + W),{2p' + W),' 



pq 



{.P + P')u 



(4.17) 









-9,iu + 



p'J2p' + 2q + W), 



p'q 





{P + P')u- 



(4.18) 



In table ^]T] the corresponding branching ratios are presented for the four decay 



modes for E^ > SOMeV and 61-y > 20°. For K^^^, the rates are to be understood 
as T{Kl -^ TT^Puj). The number of events correspond to the design values for 
DA<I>NE (cf. App. 0). 
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Table 4.1: Branching ratios for tree level amplitudes for E^ > 30MeV and 6i-y > 20 
in the K rest frame. 



decay 



-"■687 

-"■687 

^/^37 



BR(tree) 

2.8 X 10"^ 

1.9 X 10-5 
3.6 X 10-3 
5.2 X 10"^ 



#events/yr 



2.5 X 10^ 
1.7 X 10^ 
4.0 X 10^ 
5.7 X 10^ 



Table 4.2: Experimental results for Ki-^^ decays 



decay 


exp. 


-t-''y,min 


7^ events 


BR 






m 


lOMeV 


192 


(2.7 ±0.2) X 10-^ 


0.6 < cos ^^7 <0.9 


Kt,, 


47 


lOMey 


13 


(3.7 ±1.3) X 10-^ 


77 


Kt,. 


4^ 


SOMeV 


16 


(2.3 ±1.0) X 10-^ 


cos6'e^ < 0.9 


-"■^37 


45 


SOMeV 





< 6.1 X 10-5 


90% c.l. 


^°87 


m 


IbMeV 


10 


(1.3 ±0.8) X 10-2 





4.3 Previous experiments 

The data sample for Ki^^ decays is very limited and it is obvious that DA<I>NE 
will be able to make significant improvements. The present experimental status is 



summarized in table 4.2. 



A comparison between tables ^]T| and ^4.2| shows the tremendous improvement in 
statistics to be expected at DA$NE. We shall come back to the question whether 
this improvement will be sufficient to test the standard model at the next-to-leading 
order, O(p^), in CHPT. 



4.4 Theory 

Prior to CHPT, the most detailed calculations of Ki^^ amplitudes were performed 
by Fearing, Fischbach and Smith |^D[ using current algebra techniques. 

In the framework of CHPT, the amplitudes are given by ( 4.17|) and ( [4.18 ) to 
leading order in the chiral expansion. 



4.4.1 CHPT to 0(/) |n| 



There are in general three types of contributions [^ : anomaly, local contribu- 
tions due to £4 and loop amplitudes. 
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Figure 4.2: Loop diagrams (without tadpoles) for K13 at 0{p'^). For Ki^^, the photon 
must be appended on all charged lines and on all vertices. 



The anomaly contributes to the axial amplitudes 



At. = ^^^ ls,,,^q'{4p' + Wr + j^^y^^^j^e^xpaW^p'^'W 



flU 



4° 



167r2F2 
i 

~87r2F2 



W^ - Ml 



(4.19) 



Suup.q'W. 



The loop diagrams for i^^s^ are shown in Fig. [4.2| . We first write the Kj^ matrix 
element in terms of three functions fi~, f2, fit which will also appear in the invariant 
amplitudes V^^. Including the contributions from the low-energy constants L^,Lg 
in £4, the Ki^, matrix element F+ is given by 



fHt) 
ft it) 

ft it) 

L'g 

A/ 



liMl-M'^-t)fnt) + ftit) 



ftit)p'.+ 

AT -^ 

V2 + ^^t + 2Y.ici-c{)Bi{t) 



ip-p)i 



v^F2 t ^ 



1=1 
) t J J\ 



1=1 

3 



2Bi{t) 



(t + Aj)AjJi{t) 
2t 






ciAiJj{t) \ 
(4.20) 



Lli^^) 

M] -mj , t 



2567r2 



In 



iP-pT 
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Table 4.3: Coefficients for the K^ loop amplitudes corresponding to the diagrams 



/ = 1, 2, 3 in Fig. O- All coefficients cf must be divided by 6 v^F^ 



/ 


Mi 


mi 




4 


ci 


1 

2 
3 


Mk 
Mk 


M^ 

Mr, 

Mk 


1 
3 



-2 
-6 
-6 


-Ml - 2MI 
-Ml - 2Ml 

-ml 



Lq is a scale independent couphng constant and we have traded the tadpole con- 
tribution together with L5 for Fr-/F^ in fjKt). The sum over I corresponds to the 
three loop diagrams of Fig. ^^ with coefficients c{, Cj, c| displayed in table |0 . 
We use the Gell-Mann-Okubo mass formula throughout to express M^ in terms of 
Ml, Ml. The functions J/(t) and Bi{t) can be found in App. @. 



The standard Ki-^ form factors f-\-{t), /-(t) as given in the previous subsection |24 



are 



hit) 
f-{t) 



T2^^^'^ 



(4.21) 



v^ 



{Mi-Mi-t)fnt)+vt{t)~fnt) 



It remains to calculate the infrared finite tensor amplitude ICt- The invariant 
amplitudes V^ can be expressed in terms of the previously defined functions f^ and 
of additional amplitudes /i,/2,/3. Diagrammatically, the latter amplitudes arise 



from those diagrams in Fig. ^^ where the photon is not appended on the incoming 
K^ (non-Bremsstrahlung diagrams). The final expressions are 



w„ 



h + pWjtiwD + ftiw, 



p'wjt{w^,)+n{w, 



I2- — 

pq 
pq 



p'WftiW^) + UiW) - p'WJ^iW') - f^iW')] (4.22 



W + q 



pq 



p — p 
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The amplitudes /i, /2, h in Eq.( ^4.22| ) are given by 
h - 



'^^:(I^ + M + ^X-. 



V2F2' '"' V2F2 

p'W,AW^-Aj)Jj 



7=1 



(W2_A,)J, 



2W^' 



2Gj) 



Co - Ci 



1^2 
9 



iy2 

9 



+ 4S20+p'(W^-g)L^ 



^^4tI7^ [p'^(^^ - (^' + ^^)^^) + P'^(^2 - 52^)] } 



/2 
h 

Fi 



qW 

^^" +47hci-c{)\Fj-{W' + Aj)G, 



v/2F2 qW fr{ 



4Lq 



V2F2 ,tt 



+ E 2(4 - c() 






— c 



1 W^ 



(4.23) 






2567r2 
327r2A7 ""mi 



5 - 



Ml 



4 



giy 



^ 



W^Bi - W^Bi 



{W' + Aj)Jj - (W^ + Aj)Jj 



QAtt' 



Gi = ^-fC{WlW\Mlm]) + ^-^ 

Jj = Jj{W'), Ji = Jj{W;) 
Bi ^ BiiW'), Bi^B'.iW'). 



The function C{W^,W'^,Mf,mj) is given in App. ^. All the invariant ampli- 
tudes Vj^, . . . , V4~ are real in the physical region. Of course, the same is true for the 
Ki3 matrix element F+. 

The K^2 amplitude has a very similar structure. Both the K^^ matrix element F^ 
and the infrared finite vector amplitude V^^, can be obtained from the corresponding 
quantities F+ and V^ by the following steps: 

• interchange p' and —p ; 

• replace ^ by — ^ in /3+; 



insert the appropriate coefficients cj for Kf^^ listed in table |44 
multiply F+ and V^ by a factor —\/2. 
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Table 4.4: Coefficients for the Kf.^^ loop amplitudes corresponding to the diagrams 



/ = 1, 2, 3 in Fig. O. All coefficients cf must be divided by 6^2^^ 



/ 


Mi 


mi 


o[ 


4 


•-a 


1 

2 
3 


Mk 
Mk 

M^ 


Mk 



6 

4 


-3 
-3 

-2 


-3M| 

Mi + 2M2 
-2Mi + 2Ml 



Table 4.5: Branching ratios and expected number of events at DA$NE for K^ 



VS-y- 



^e37 


BR 


7^events/yr 


full O(p^) amplitude 

tree level 

O(p^) without loops 


3.0 X 10"^ 
2.8 X 10"^ 
3.2 X 10"^ 


2.7 X 10^ 
2.5 X 10^ 
2.9 X 10^ 



K.. 


BR 


#events/yr 


full 0{p^) amplitude 

tree level 

0{p^) without loops 


2.0 X 10~5 
1.9 X 10-5 

2.1 X 10-5 


1.8 X 105 

1.7 X 105 

1.9 X 105 



4.4.2 Numerical results 



In calculating the rates with the complete amphtudes of the previous subsection, 



we use the same cuts as for the tree level rates in Subsect. 4.2 



E^ > 30MeV 
Oi^ > 20°. 



(4.24) 



The physical values of M^ and M^ are used in the amplitudes. M^ is calculated 
from the Gell-Mann-Okubo mass formula. The values of the other parameters can 
be found in Ref. and in appendix ^ 

The results for Kj^^ and K^^^ are displayed in tables ^^ and ^^, respectively. 
For comparison, the tree level branching ratios of table ^]l] and the rates for the 
amplitudes without the loop contributions are also shown. The separation between 
loop and counterterm contributions is of course scale dependent. This scale de- 
pendence is absorbed in the scale invariant constants Lg,LiQ defined in Eqs.( [4.20l) , 



( [4.23| ). In other words, the entries in tables [4.5| , ^.6| for the amplitudes without loops 
correspond to setting all coefficients cf in tables ^^, ^^ equal to zero. 
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Table 4.6: Branching ratios and expected number of events at DA$NE for i^fg 



«37- 



^e37 


BR 


#events/yr 


full 0{p'^) amplitude 

tree level 

0(j9*^) without loops 


3.8 X 10-3 
3.6 X 10-3 
4.0 X 10-3 


4.2 X 10^ 
4.0 X 10^ 
4.4 X 10^ 




^m37 


BR 


#events/yr 


full O(p^) amplitude 

tree level 

0(p'^) without loops 


5.6 X 10-^ 
5.2 X 10-^ 
5.9 X 10-^ 


6.1 X 10^ 

5.7 X 10^ 
6.5 X 10^ 



4.5 Improvements at DA$NE 

The numerical results given above demonstrate very clearly that the non-trivial 
CHPT effects of O(p^) can be detected at DA<I>NE in all four channels without any 
problem of statistics. Of course, the rates are bigger for the electronic modes. On 
the other hand, the relative size of the structure dependent terms is somewhat bigger 
in the muonic channels (around 8% for the chosen cuts). We observe that there is 
negative interference between the loop and counterterm amplitudes. 

The sensitivity to the counterterm coupling constants Lg, Liq and to the chiral 
anomaly can be expressed as the difference in the number of events between the tree 
level and the O(p^) amplitudes (without loops). In the optimal case of K^^^, this 
amounts to more than 4 x 10^ events/yr at DA$NE. Almost all of this difference is 
due to Lg. It will be very difficult to extract the coupling constant Lio from the total 
rates. A more detailed study is needed to determine whether Lio can be extracted 
from differential distributions. 



The chiral anomaly is more important for K^ , but even there it influences 
the total rates rather little. Once again, a dedicated study of differential rates is 
necessary to locate the chiral anomaly, if possible at all. 

On the other hand, taking into account that Lg is already known to good accuracy 
(see table 1 in Ref. 0), Ki^^ decays will certainly allow for precise and unambiguous 
tests of the one-loop effects in CHPT |TT|| . 
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5 Ki4 decays 

In this subsection we discuss the decays 

K+(p) -^ n+ipi)n-ip2)l+ipi)uiipu) (5.1) 

K+{p) ^ 7r°(pi) 7r°(p2) /+(P/) ///(J9.) (5.2) 

K\p) ^ n\p,)n-{p,)l+{pi)ui{p,) (5.3) 

and their charge conjugate modes. The letter I stands for e or fi. We do not consider 
isospin violating contributions and correspondingly set m^ = rrid, a = 0. 

5.1 Kinematics 

We start with the process (|5.1j ). The full kinematics of this decay requires five 
variables. We will use the ones introduced by Cabibbo and Maksymowicz [^. It is 
convenient to consider three reference frames, namely the K^ rest system (S/^), the 
TT+TT^ center-of-mass system (S27r) and the l^ui center-of-mass system (S^j^). Then 
the variables are 

1. Stt, the effective mass squared of the dipion system, 

2. si, the effective mass squared of the dilepton system, 

3. 6t^, the angle of the vr"*" in S2,r with respect to the dipion line of flight in Hk, 

4. 9i, the angle of the /+ in S^,^ with respect to the dilepton line of flight in T,^, 
and 

5. (f), the angle between the plane formed by the pions in J^x and the correspond- 
ing plane formed by the dileptons. 



The angles 9^^, 9i and (p are displayed in Fig. |5.1| . In order to specify these 
variables more precisely, let pi be the three-momentum of the tt"'" in J]2n and pi the 
three- momentum of the l~^ in S;,^. Furthermore, let w be a unit vector along the 
direction of flight of the dipion in S/^, and c{d) a unit vector along the projection 
oi pi{pi) perpendicular to v{—v), 

c = (pi - vv ■ Pi)/[Pi - (pi ■ tJ) Y^^ 
d = {pi-vv- pi)/[pi^ - {pi ■ vfY^"^ . 



The vectors w, c and d are indicated in Fig. |5.1| . Then, one has 

cos9^ = v-pi/\pi\, cos9i = ~v-pi/ \pi\ 
COS0 = c-d, sinKJ) = {ex v) ■ d. (5.4) 
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Figure 5.1: Kinematic variables for i^^ decays. The angle 9.,^ is defined in E27r, di in 

Sjj, and in Sx- 

The range of the variables is 

4M2 < s^< {Mk - mif 



m] < si< {Mk - v^)' 
< 9^,01 < 71,0 <(/)< 27V. (5.5) 

It is useful to furthermore introduce the following combinations of four vectors 



P = Pi+P2, Q = Pi-P2 , L = pi+p^ , N = pi-p^ 
together with the corresponding Lorentz invariant scalar products 

p2 = g^^ Q^ = AMl-s^, L^ = si, N^ = 2mf-si, 
PQ = 0, 



PL 

PN 
QL 
QN 

LN 
< LNPQ > 



-(Ml-s^-si), 
ziPL + {l - zi)Xcosei, 

(T^XCOS^TT) 

ZiQL + (T7r(l — Zl) [P L COsOt, COS 9l 

— (sttSiY'^ sin 9t, sin 9i cos > 



m, 



-fiupa 



L'^N'^PfQ'' 

a/2 



— (^St^Si) ' aj^(l — Zi)X sin 9j^ sin 9i sin I 



with 



X = {{PLf-s^siY/' = h'/\Ml,s^,si) 



Zl 



(y-n 



mf/si 



(5.6) 



(5.7) 



(5.^ 
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Below we will also use the variables 

t = (pi - pf 

u = {p2-pf. (5.9) 

These are related to St,, si and 9t, by 

t + u = 2Ml + Ml + si-s^ 

t-u = -2a^Xcos9^ . (5.10) 

5.2 Matrix elements 

The matrix element for K~^ —>■ -k^-k^I'^vi is 



G 



F 



where 



T = -^V:MPuh,il - iMpiW^ - An (5.11) 



I^. = < n+{p,)n-{p2)ont \ I^^^'^O) \ K+{p) >; I = V,A 

TT 

A^ = -i-^[P^F + Q^G + L^R] (5.12) 



and eoi23 = 1- The matrix elements for the other channels (|5.2|J5.3| ) may be obtained 
from ( |5.11|j5.1^ ) by isospin symmetry, see below. 



The form factors F, G, R and H are real analytic functions of the three variables 
P1P2, PiP and p2p. Below, we will use instead the variables {s-,,, si, 9,,} or {s-,,, t, u}. 

Remark: In order to agree with the notation used by Pais and Treiman [^ and 
by Rosselet et al. |^, we have changed our previous convention [^, ^ in the 
definition of the anomaly form factor H. See also the comments after Eq. ( |5.21| ). 



5.3 Decay rates 

The partial decay rate for ( |5.1|) is given by 



c^r = ^^ .^ ,g II I ^ P dLiPs{p;Pi,Pu,Pi,P2)- (5.13) 



spins 



The quantity Y^gpi^s | T p is a Lorentz invariant quadratic form in F,G,R and 
H. All scalar products can be expressed in the 5 independent variables St,,si,9.„, 9i 
and (/), such that 

J2 \T\^= 2G| I Vus l^M^^Ms^, si, 9^, 9i, 0) . (5.14) 



spins 
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Carrying out the integrations over the remaining 4-3 — 5 = 7 variables in ( |5.13| ) 
gives 1^ 



dT^ = Gp I Vus 1^ N{s^,si)J^{s^,si,9^,9u(j))ds^dsid{cos9^)d{cos8i)d(f) (5.15) 

where 

N{s^, si) = (1 - zi)a^X/{2^^T:''Ml) . (5.16) 

The form factors F, G, R and H are independent of and 6i. It is therefore possible 
to carry out two more integrations in ( |5.15|) with the result 



dV^ = G% I Vus 1^ N{s^,si)J3{s^,si,9^)ds^dsid{cos9^ 
The explicit form of J5 is 



(5.17) 



J. 



{PLy-(PNf-s^s, + mfs^ 



2^2 



+ |G|' {QLy-iQNy-Q'si + mfQ 



|2 _2 



1 



+ 



Mi 



ifl 



Si -nil 
{m^ - Si) Q^X^ + s^{QLf - < LNPQ >^ 



+ {F*G + FG*) [{PL) (QL) - (PN) (QN)] 
+ {F*R + FR*) mf [{PL) - (PN)] 

+ -T7TiP*H + FH*) \{QN){PLf-{QL){PL){PN)-sMQN)+mls^{QL) 



M' 



K 



+ {G*R + GR*)mi[{QL)-{QN)] 



+ -^{G*H + GH*) UpL){QL){QN) - {PN)iQL)^ + si{PN)Q^ - m^iPL)Q' 
+ TPT < LNPQ > \-{F*G - FG*)Ml + {F*H - FH*){PN) 



Ml 



{G*H - GH*){QN) + {R*H - RH*)mi 



(5.18) 



For data analysis it is useful to represent this result in a still different form which 
displays the 9i and (p dependence more clearly ||5^ : 



J5 = 2(1 — zi) Ii + I2 cos 29i + J3 sin^ 9i ■ cos 20 + I4 sin 29i ■ cos + /s sin 9i ■ cos < 
+ Jfi cos ^i + I7 sin 6'i ■ sin d) + Ir sin 26'/ ■ sin (p + Iq sin^ ^; ■ sin 2 



(5.19) 



One obtains 



/i = \[{l + zi)[F,[' + ^{3 + zi){[F,[' + [Fs[')sm'9^ + 2zi[F,['^ 
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2 — 


-\il-zi)[\F,^-l{\F,^ + \Fs')sin^e 


3 = 


4 '^ J 


4 = 


^il-zi)Re{F*F2)sme^ 



h 
h 
h 
h 



- { Re(F;F3) + zi Re{FlF2)} smO^ 

- { Re(F2*F3) sin^ 9^ - zi Re{F*F^)] 

- { ImiF^F^) + zi lm{F:F,)}sme^ 

^{1- Zi)lm{F*Fs) sine^ 
-l-{l-Zi)lm{F;F,)smH^, 



(5.20) 



where 



Fi 

F2 
F3 

F, 



X-F + a^{PL)cose^-G 



Ml 



-{PL)F - siR - a^X cos 6^ ■ G 



(5.21) 



The definition of Fi, . . . , F4 in (|5.21|) corresponds to the combinations used by Pais 
and Treinian [^ (the different sign in the terms ~ PL is due to our use of the 

metric diag(H )). The form factors /i, . . . , /g agree with the expressions given 

in |5^. We conclude that our convention for the relative phase in the definition of 
the form factors in Eq. (|5.12|) agrees with the one used by Pais and Treiman. The 
comparison of ( p. IS] ) with [^, table II] shows furthermore that it also agrees with 
this reference. 



The quantity J3 can now easily be obtained from ( |5.19| ) by integrating over (j) 
and 61, 

J3 = [ d<i) d{cosei)J5 = 87r(l - zi) III - -I2] . (5.22) 



5.4 Isospin decomposition 



The Kii decays (|5.2|) and ( |5.3| ) involve the same form factors as displayed in Eq. 
( ^.12|) . We denote by A+_, Aqq and Aq_ the current matrix elements of the processes 
(|5.1|)-(p73D. These are related by isospin symmetry 0, 



A, 



Ao- 

V2 



A 



00 



(5.23) 



'We use the Condon-Shortley phase conventions. 
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This relation also holds for the individual form factors, which may be decomposed 
into a symmetric and an antisymmetric part under t ^-* u {pi ^-* p2). Because of 
Bose symmetry and of the A/ = | rule of the relevant weak currents, one has 

{F,G,R,H)oo = —{F^,G~,R^,H )+_ 

{F,G,R,H)o_ = V2iF-,G+,R-,H+)+_ (5.24) 

where 

F±_ = ^ [F{s^, t, u) ± F(s^, M, t)] (5.25) 

and F{sT,,t,u) is defined in Eq. ( p.l2| ). 

The isospin relation for the decay rates is 

r{K+ -* n+7i-l+ui) = ^^{Kl -^ 7r\^l^u) + 2r{K+ -^ n\^l+ui) . (5.26) 

Isospin violating contributions affect the matrix elements and phase space, as a result 
of which this relation is modified. In order to illustrate the (substantial) effects from 
asymmetries in phase space, we take constant form factors F, G and set R = 0,H = 
0. Eq. ( p. 26 ) then reads (with physical masses for K'^ -^ 7r"'"7r~/"'"z/;,7r°7r°/"'"i/i and 
with M^o = M^± = 137 MeV in Kl -^ 7r°7rT/±z/) 

(16.0^2 + 3.1^^)10 = (20.1^2 + 2.0G2)ro 

To = KVlO'sec-^ (5.27) 

in the electron mode and 

(1.79F2 + 0.25G'2)ro = (2.64^2 + 0.20^^)10 (5.28) 

in the muon mode. 

5.5 Partial wave expansion 

The form factors may be written in a partial wave expansion in the variable 6.,^. We 



consider a definite isospin vrvr state. Suppressing isospin indices, one has ||56|, ^ 



F = Y,Pi{cose^)fi-^^cose^G 

1=0 ^ 



G = ^P/(cose.)(7; 



1=1 
R = '^ Pi{cos 6-^)ri + ^^ cos 6-^G 

1=0 ^ 

oo 

H = ^P/(cos^^)/ii (5.29) 

1=0 
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Table 5.1: Rates of K^ decays [||. The data for Kl -^ Tr^vr^e^z/ is from ref. |^ 





tl events 






branching 
ratio 


Particle Data 
Group 


DA<1>NE 

1 yr 


improve- 
ment 


K^ —^ TT + TT e+Z/e 

K+ -^ TrOvrOe+z/e 

Kl -^ 7r°7r±e^z/ 

K~^ —>■ 7r+7r-yU"^z/^ 


3.91-10-5 

2.1-10-5 

5.16- 10-5 

10-6 


3-10^ 
<50 

729 

7 


3-105 
2-105 
7-10^ 
2-10^ 


10 

> 4 - 10^ 

7-102 

3 -103 



where 



/>;(.) = |p,(.) . 



(5.30) 



The partial wave amphtudes fi,gi,ri and hi depend on s^ and s;. Their phase 
coincides with the phase shifts Sj in elastic tttt scattering (angular momentum I, 
isospin /). More precisely, the quantities 



e-''''X2i 



e-^^^mX2m; / = 0,1,... ; X = /,^,r,/i 



(5.3i; 



are real in the physical region of K14 decay. The form factors Fi and F4 therefore 
have a simple expansion 



Fi = xY,Pi{cose^)fi 
I 

^4 = -J2Pi{cose^){PLfi + siri). 

I 



(5.32) 



On the other hand, the phase of the projected amphtudes 
F21 = I Piicos9^)F2dicose^) ; / = 0, 1, 2, . . . 



(5.33) 



is not given by (5/, e.g., e *'^i-F2o is not real in the isospin one case. A similar remark 
applies to F3. 



5.6 Previous experiments 



We display in table p.l| the number of events collected so far. The data are obviously 
dominated by the work of Rosselet et al. [Q, which measures the tt+tt- final state 
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with good statistics. The authors parametrize the form factors as 

F = f,e< + fpe'^'^ cos e^ + D-wme 

G = ge'^^^ + D-wave 

H = /le'"^! + D-wave (5.34) 

with fs,fp,g and h assumed to be real f\ . Furthermore, they put nie = 0, such 
that the form factors R and F4 drop out in the decay distribution. Despite the 
good statistics, the experiment has not been able to separate out the full kinematic 
behaviour of the matrix elements. Therefore certain approximations/assumptions 
had to be made. For example, no dependence on si was seen within the limits of the 
data, so that the results were quoted assuming that such a dependence is absent. 
Similarly, fp was found to be compatible with zero, and hence put equal to zero 
when the final result for g was derived. A dependence on Stt was seen, and found to 
be compatible with 



fs{q') = 


/,(0)[1 + A;g2] 


9{<f) = 


g{0)[l + Xgq'] 


Hq') = 


h{0)[l + Xhq^] 


q' = 


{s^- AMD /AMI 


A/ = 


Xg = Xh = X. 



(5.35) 

with 

(5.36) 

These approximations to the form factors do not agree completely with what 
is found in the theoretical predictions. Dependence on si and non-zero values for 
higher partial waves all occur in the theoretical results. 

The experimental results for the threshold values and the slopes of the form 
factors are |^ 

/,(0) = 5.59 ±0.14 

^(0) = 4.77 ±0.27 

h{0) = -2.68 ±0.68 

A = 0.08 ±0.02. (5.37) 

We have used | Vus \= 0.22 in transcribing these results. (We note that from Eqs. 
(p:^-|07|) and fp = Owe obtain Tk,^ = (2.94±0.16)-103 sec'^ This value must be 
compared with T^^^ = (3.26 ±0.15) • 10^ sec~^ obtained in the same experiment.) In 
addition to the threshold values (|5.371 ) of the form factors, the phase shift difference 
6 = 6q — 6l was determined |53| in five energy bins. The S-wave scattering length 



Oq was then extracted by using a model of Basdevant, Froggatt and Petersen |]58 



6 



Note that, according to what is said in the previous subsection, the terms denoted by "D-wave" 



in Eq. ( |5.34| ) aU contain (complex) contributions which are proportional to Pj(cos0^),/ > 0. 
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This model is based on solutions to Roy equations. The result for the scattering 
length is 

(5.38) 



0.28 ±0.05. 



A study by [|^, based on a more recent solution to Roy equations, gives 



0.26 ±0.05. 



(5.39) 



Turning now to the Tr^Tr^e^z/g channel, we consider the following recent data 
(based on 729 ± 15 events) [|60|] : 



BR{Kl -^ TT^TT^e^u) = (5.16 ± 0.2 ± 0.22) ■ 10"^ 
The group also measured the G form factor. Defining 



Go^ = Go(l ± A,g2)e*^i 



(5.40) 



(5.41) 



they find 



Go 
A„ 



7.8 ±0.7 ±0.2 , 
0.014 ±0.087 ±0.070 . 



(5.42) 



The slope agrees within the errors with the value (|5.37|) found by Rosselet et al. ||53| . 
To compare the value of the form factor at threshold, we use the isospin prediction 



|Go|/V2 = 5.5±0.5 



(5.43) 



which is not incompatible with g{0) = 4.77 ± 0.27 in eq. (|5.37| ). (Here we have used 
\Vus\ = 0.22 to transcribe the data. Furthermore, we assume that the form factor 
Go- measured in Ref. ||60| indeed has to be divided by -\/2 for the comparison with 
|]53[] . This is not quite clear to us reading [pO| .) 

Finally for the channel TT^TT^e'^Ue, we consider the rate P] 



K+^TrOirOe+u, 



(l.70+°;3^) . lO^sec-^ 



(5.44) 



The kinematic dependence of the form factors on the variables s^, si and 6.^ has 
not yet been resolved experimentally in this decay. In order to proceed, we assume 
that the form factors in this channel are independent of 9.„, e.g., Foo = ^00(5,^,^ ± 
u) etc. As a result of this assumption. Goo and Hqq vanish by Bose statistics. 
The contribution from Rqq is completely negligible in the electron mode, and the 
contribution from the anomaly form factor to the decay ( 15.44] ) is tiny. We neglect 
it altogether, as a result of which the above decay is fully determined by Foo. We 
write 

Foo = Fo{l + Xq^)e< , (5.45) 

and obtain for the rate 
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K+^TvOTvOe+iye 



|FoV;.r(2.01 ± 1.7A ± 0{X^)) ■ lO^sec^^ . 



(5.46) 



51 



This finally gives with A = 0.08 



\Fn\ = 5.72" 



rO.57 



0| — ^. 1^-0.49 5 

which compares very well with the isospin prediction 

|i"o| = |/s(0)| = 5.59 ±0.14 



(5.47) 



(5.48) 



5.7 Theory 

The theoretical predictions of K^ form factors have a long history which started 
in the sixties with the current algebra evaluation of F, G, R and H. For an early 
review of the subject and for references to work prior to CHPT we refer the reader 
to [^ (see also [^). Here we concentrate on the evaluation of the form factors in 
the framework of CHPT 0, Em Fl, B2]. 



5.7.1 Form factors at tree level 



The chiral representation of the form factors at leading order was originally given 
by Weinberg [|6^] , 



G 



M 



K 



R 
H 



V2F^ 



2V^F, \si-Ml 
. 



3.74 , 
+ 1 



(5.49) 



The next-to-leading order corrections are displayed below, and the later sections 
contain an estimate of yet higher-order contributions. Here we note that the total 
decay rates which follow from Eq. ( [5.49| ) are typically a factor of two (ore more) 
below the data. As an example, consider the channel K^ — >■ 7r^7r~e^t'e- Using 
( [5.49[ ), the total decay rate becomesQ 1297 sec~^, whereas the experimental value is 
3160±140 sec-i g. 



5.7.2 Form factors at one loop 



The one-loop result for F [^,||55| may be written in the form 
F(s^, t, u) = -^ Ii + 1-(Uf + Pp + Gf)+ 0{E') 



(5.50) 



■^If not stated otherwise, we use F^ = 93.2 MeV, |Ks| = 0.22 and {AU,Mk) = (139.6,493.6) 
MeV, (135,493.6) MeV and (137,497.7) MeV for the decays (|5^), (fj) and (fj), respectively. 
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The contribution UF{sn,t,u) denotes the unitarity correction generated by the one- 
loop graphs which appear at order E'^ in the low-energy expansion. Its expression 
will be given in appendix 0. 

The contribution Pp{s.,^,t,u) is a polynomial in STj,t,u obtained from the tree 
graphs at order E^. We find 



where 



PF{s-„,t,u) =^Pi^F{s-K,t,u)Ul , 



i=l 



,F 


= 


32(s^-2M2) , 




,F 


= 


8{Ml + s^-si) , 




,F 


= 


A{Ml - 3M^ + 2s^ 


-t) 


,F 


= 


32M2 , 




,F 


= 


4M,2 , 




.F 


= 


2si . 





(5.51) 



(5.52) 



The remaining coefficients Pi^p are zero. 

The contributions Cp contain logarithmic terms, independent of s.^,t and u: 



C, 



1 



2567r2 



/l/f2 /l#2 /1#2 

5M^Mn^-2M|,ln^-3Mjln^ 



(5.53) 



The corresponding decomposition of the form factor G is ||54||,||55 

G{s^,t,u) 



^"^ h + l-(UG + PG + CG)+0{E') 



V2E^ I F^ 



For the expression of Ug see appendix 0. The polynomials 



Pg = ^Pt^G{.s^,t,u)Ll 



i=l 



are 



P2,G = 8{t-u) , 

P3,G = A{t-Ml-Ml) 

p,,G = AMI , 

P9,G = 2Si , 



(5.54) 



(5.55) 



(5.56) 
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The remaining pj g vanish. The logarithms contained in Cq are 

Cg = —Cp- 



(5.57) 



The form factor R contains a pole part 2'(s^, t, u)/{si — M"^) and a regular piece 
Q. [Since the axial current acts as an interpolating field for a kaon, the residue of 
the pole part is related to the KK -^ vrvr amplitude in the standard manner.] We 
write 



R 
I 



M 



K 



2V2F^ U, - M 



2 rQ + OiE^) 

K 



Bj + —{Uj + Pj + Ci), I = Z,Q 



(5.5J 



According to ( ^.491) , the Born terms Bj are |6 



Bz 
Bo 



u 



(5.59) 



The one-loop corrections have been worked out in Ref. [£2|. The unitarity correc- 
tions Uj are displayed in the appendix 0. The residues Pz and Cz are 



with 



Pz{s^,t,u) = J2Pi,zis^,t,u)P 



r 
i 1 



i=l 



(5.60) 



,z 


= 


■62[ 


s^-2M^)[s^- 


2m;) , 




,z 


= 


8{sl + U-) , 




,z 


= 


-2 


'2{i^ + AJ:)s^-5sl-u^- 


- IQMlMi 


,z 


= 


32 


'T.s^-AMlMl 


; 




,z 


= 


4[( 


s^ + z/)s - mlMl 


7 




,z 


= 


i2mlMl , 




,z 


= 


64i 


Km'. , 









E = Ml + M^ . 
The remaining pi^z vanish. The logarithms in Cz are 



c. 



Ml 



1287r2 
For the nonpole part Q we find: 



3M^ln^-2Miln^-M,5lnqL 



(5.61) 



(5.62) 
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PQ{s^,t,u) = ^pi^Q{s^,t,u)Ll 



(5.63) 



with 



Pi,Q = 32(s^-2M^) 



P2,Q 
P3,Q 
P4,Q 
P5,Q 
P9,Q 



4(s. - 2M2) + M| - s. 



32Mi 



4S 



(s^ + z/) - (M^ - si) 



The remaining p, q vanish. The logarithms in Cq are 

1 



Cr 



1287r2 



A/f2 /\/f2 M' 

5M,Mn ^ - 2Mi In ^ - 3M^Mn ^ 
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(5.64) 



(5.65) 



The first nonvanishing contribution in the chiral expansion of the form factor H 
is due to the chiral anomaly IT^. The prediction is 



H 



V2MI 



K 



87r2F3 



-2.66 



in excellent agreement with the experimental value |5^ H = —2.68 ± 0.68 
next-to-leading order corrections to H have also been calculated |]65 



(5.66) 

The 
If the new 



low-energy parameters are estimated using the vector mesons only, these corrections 
are small. 

The results for F, G and R must satisfy two nontrivial constraints: i) Unitarity 
requires that F, G and R contain, in the physical region 4M^ < s^^ < {Mk — rniY, 
imaginary parts governed by S- and P-wave vrvr scattering [these imaginary parts 
are contained in the functions Ao(s^), Ai(s^)]. ii) The scale dependence of the low- 
energy constants L[ must be compensated by the scale dependence of Uf,g,z,q and 
Cf,g,z,q for all values of St^^I^u^M'^^M'^. [Since we work at order E'^ in the chiral 
expansion, the meson masses appearing in the above expressions satisfy the Gell- 
Mann-Okubo mass formula.] We have checked that these constraints are satisfied. 

Because the one-loop contributions are rather large, one expects still substantial 
corrections from higher orders. In the following section, we therefore first estimate 
the effects from higher orders in the chiral expansion, using then this improved 
representation for the form factors in a comparison with the data. 
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5.7.3 Form factors beyond one loop 

To investigate the importance of higher-order terms, we employ the method 
developed in Ref. [Q. It amounts to writing a dispersive representation of the 
partial wave amplitudes, fixing the subtraction constants using chiral perturbation 
theory. Here, we estimate the higher-order terms in the 5'-wave projection of the 
amplitude Fi, 

fis^,si) = (AnXy' J dnF,{s^,t,si) , (5.67) 

because this form factor plays a decisive role in the determination of L^,L'2 and 
L3, and it is infiuenced by S'-wave vrvr scattering which is known |6^ to produce 
substantial corrections. 



5.7.4 Analytic properties of partial waves 

Only the crossing-even part 

F+ = XF+ + a^{PL) cos 6^ ■ G~ (5.6J 



contributes in the projection ( |5.67| ). The partial wave / has the following analytic 
properties: 

1. At fixed s;, it is analytic in the complex s^r-plane, cut along the real axis for 
Re s^ > AMI and Re s^ < 0. 

2. In the interval < s^ < 4M^, it is real. 

3. In 4M^ < Stt < 16M^, its phase coincides with the isospin zero S'-wave phase 
(5q in elastic vrvr scattering, 

U = e^<f- J± = fis^±te,si). (5.69) 



The proof of these properties is standard |^8[. Here we only note that the presence 



of the cut for s-j^ < follows from the relations 

t = M^ + ^^ ^^' ~ '- - a.Acos^. , 
t(cos^^ = -1, s^ < 0) > {Mk + M^f. (5.70) 

Since F^ and G~ have cuts at t > (M^ + M^r)^, the claim is proven. 

5.7.5 Unitarization 
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We introduce the Oinnes function 



r2(s^) = exp 



ds 6^{s) 



TC JiM'^ S S — St, 



(5.71) 



where A will be chosen of the order of 1 GeV below. According to ( |5.69D , multipli- 
cation by VL~^ removes the cut in / for 4M^ < 3,^ < 16M^. Consider now 

f = fL + fR, (5.72) 

where fLifR) has only the left-hand (right-hand) cut, and introduce 

v = n-\f-h). (5.73) 

Then v has only a right-hand cut, and we may represent it in a dispersive way. 



Vo + ViStt H 



dsimn-^if-h) 



TT J4Af2 S^ 



S - Sr 



(5.74) 



We expect the contributions from the integral beyond 1 GeV^ to be small. Further- 
more, fi~^/ is approximately real between 16M^ and 1 GeV^, as a result of which 
one has 



^^0 + ViStt 



sl f^' dsfLlmn-^ 



71 J4M2 S^ 



For given vq, f i, /l and Q, the form factor / is finally obtained from 



(5.75) 



(5.76) 



The behaviour of Jl at Stt — ^ is governed by the large- |t| behaviour of F~^ and 
G~, see ( p.7U| ). Therefore, instead of using CHPT to model fi, we approximate the 
left-hand cut by resonance exchange. To pin down the subtraction constants vq and 
Vi , we require that the threshold expansion of / and /chpt agree up to and including 
terms of order 0{E'^). For a specific choice of /i, this fixes Vq,Vi in terms of the 
quantities which occur in the one- loop representation of F^ and G~ . In the case 
where fL = 0,f has then a particularly simple form at si = 0, 



/(s^,s/ = 0)1/^=0 

^^0 



fi(f + ViStt) , 

Mk r 1 

^^1.05 + — 
V2FA F^ 



2 rr 



- QAM;L\ + 8M^L^ 



-2(Mi - ml)L, + -{Ml - AMl){AU, + L,) 



Vl 



Mk r 1 

^^0.38 + — 

V2FA F^ 



2 Ml 



32LI + 8L^ + IOL3 

■AM? 



AMI 



-{AL\ + U 



(5.77) 
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The details of the evaluation of fi, vq and of vi may be found in Ref. |]62 



In the partial wave /, the effects of the final-state interactions are substantial, 
because they are related to the / = 0, S'-wave vrvr phase shift. On the other hand, 
for the leading partial wave in G~^ = ge^^^ + ■ ■ ■, these effects are reduced, because 
the phase 5p is small at low energies. We find it more difficult to assess an estimate 
for the higher-order corrections in this case - we come back to this point in the 
following sections. 

We add a remark concerning the choice of the subtraction point in the Omnes 
function ( ^.71| ). To remove the right-hand cut in /, the modified Omnes function 



^(stt, So) = exp ■ 



So) 



TT 



ds 5n{s 



o\ 



4Af| S — So S — STT 



would do as well, with any value of Sq < 4M^. To illustrate the meaning of So, we 
consider for simplicity the form factor / in the case where /l = 0. The choice So = 
4M^ then amounts to the statement that, at threshold, there are no contributions 
from two loops and beyond by fiat. We consider this to be rather unlikely. (For a 
different opinion see ||B^.) On the other hand, we have checked that our results do 
not vary significantly if sq is taken of the order of the pion mass squared, |so| < M^. 



5.8 Determination of Li, L2 and L3 

To illustrate the usefulness of forthcoming more accurate Ki4^ data, we determine 
here the low-energy constants L[, L^ ^^d L3 from data on K^ —* 7r"^7r~e"'"z/e decays 
and on vrvr — >■ tttt threshold parameters, using the improved S'-wave amplitude / set 
up above. For a comparison with earlier work [Q we refer the reader to Ref. ||62|| . 

We perform fits to fs{0), Xf,g{0) as given in (|5.37|) and to the tttt threshold 
parameters listed in table Ol. We introduce for this purpose the quantities 






3 r 

OTT J 



\fis^,si)\ 



(5.78) 



where the factor 3/2 sin^ 9 appears because G is expanded in derivatives of Legendre 
polynomials. Below, we identify [/s(0), ^'(0)] with [f{4M^,si), g{4M^,si)], which 
depend on s/. Furthermore, we compare the slope A/ with 



Xf{s^,si) 



f{s^,si)-f{4Mlsl) AMI 



f{^Mlsi) 



Sn - 4M2 



(5.79) 



which depends on both s.^ and si. We use these dependences to estimate systematic 
uncertainties in the determination of the low-energy couplings. [In future high- 
statistics experiments, the srdependence of the form factors will presumably be 
resolved. It will be easy to adapt the procedure to this case.] 
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Table 5.2: Results of fits with one-loop and unitarized form factors, respectively. 
The errors quoted for the L['s are statistical only. The L[ are given in units of 
10~^ at the scale /i = Mp, the scattering lengths af and the slopes bj in appropriate 
powers of M^+ . 





KeA dat 


a alone 


Kei and 


TTvr data 


experiment 




one-loop 


unitarized 


one-loop 


unitarized 





LI 


0.65 ±0.27 


0.36 ±0.26 


0.60 ±0.24 


0.37 ±0.23 




LI 


1.63 ±0.28 


1.35 ±0.27 


1.50 ±0.23 


1.35 ±0.23 




L, 


-3.4 ±1.0 


-3.4 ±1.0 


-3.3 ±0.86 


-3.5 ±0.85 




al 


0.20 


0.20 


0.20 


0.20 


0.26 ±0.05 


hi 


0.26 


0.25 


0.26 


0.25 


0.25 ±0.03 


-10 al 


0.40 


0.41 


0.40 


0.41 


0.28 ±0.12 


-10 bl 


0.67 


0.72 


0.68 


0.72 


0.82 ±0.08 


10a} 


0.36 


0.37 


0.36 


0.37 


0.38 ±0.02 


lO^fti 


0.44 


0.47 


0.43 


0.48 




lO^a" 


0.22 


0.18 


0.21 


0.18 


0.17 ±0.03 


lO^al 


0.39 


0.21 


0.37 


0.20 


0.13 ±0.3 


xVN^o. 


0/0 


0/0 


8.8/7 


4.9/7 





We have used MINUIT |£0[ to perform the fits. The results for the choice si = 
0,St, = AAMl are given in table |5^ . In the columns denoted by "one-loop", we 
have evaluated /, g and A/ from the one-loop representation given above ^ In the 
fit with the unitarized form factor (columns 3 and 5), we have evaluated / from Eqs. 
( ^.76| ), inserting in the Omnes function the parametrization of the tttt S'-wave phase 



shift proposed by Schenk |72, solution B]. For the form factor G, we have again used 
the one- loop representation. The statistical errors quoted for the Lj's are the ones 
generated by the procedure MINOS in MINUIT and correspond to an increase of 
X^ by one unit. 

A few remarks are in order at this place. 



1. It is seen that the overall description of the vrvr scattering data is better using 
the unitarized form factors, in particular so for the D-wave scattering lengths. 

2. The errors quoted do not give account of the fact that the simultaneous deter- 
mination of the three constants produces a strong correlation between them. 
To illustrate this point we note that, while the values of the Lj's in column 4 
and 5 are apparently consistent with each other within one error bar, the x^ 
in column 5 increases from 4.9 to 30.7 if the Lj's from column 4 are used in 



^We always use for L 



4i ■ 



, Lg the values quoted in table 1 of Ref. 
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the evaluation of x^ in column 5. (For a discussion about the interpretation 
of the errors see IffT 



3. The low-energy constants li, I2 which occur in SU{2)i x SU{2)ji analyses may 
be evaluated from a given set of L^, L2 and L3 |^^. Their value changes in 
a significant way by using the unitarized amplitude instead of the one-loop 
formulae: the values for (/i, I2) in column 4 and 5 are (—0.5 ± 0.88, 6.4 ± 0.44) 
and (-1.7 ± 0.85, 6.0 ± 0.4), respectively. 



4. L[, L2 and L3 are related to vrvr phase shifts through sum rules [^ Q. In 
principle, one could take these constraints into account as wel]0. We do not 
consider them here, because we find it very difficult to assess a reliable error 
for the integrals over the total vrvr cross sections which occur in those relations. 

The statistical error in the data is only one source of the uncertainty in the low- 
energy constants, which are in addition affected by the ambiguities in the estimate 
of the higher-order corrections. These systematic uncertainties have several sources: 

i) Higher-order corrections to g have not been taken into account. 

ii) The determination of the contribution from the left-hand cut is not unique. 

iii) The quantities / and g depend on si, and Xj is a function of both si and s,r- 

iv) The Omnes function depends on the elastic vrvr phase shift and on the cutoff 
A used. 



We have considered carefully these effects 62 , and found that the best determi- 



nation of L^, L2, and L3 which takes them into account is 



lO^L^M,) = 


-- 0.4 ±0.3 


103L^(M,) = 


= 1.35 ±0.3 


lO^L. = 


= -3.5 ±1.1 



(5.80) 



These values are the ones quoted in table 1 in Ref. 0. For SU{2)i x SU{2)ji 
analyses it is useful to know the corresponding values for the constants li and T2, 

h = -1.7 ±1.0 , 

k = 6.1 ±0.5 . ^^-^^^ 

The value and uncertainties in these couplings play a decisive role in a planned 
experiment [^ to measure the lifetime of tt+tt^ atoms, which will provide a com- 
pletely independent measurement of the vrvr scattering lengths laj] — OqI. 

^We thank B. Moussallam for pointing this out to us. 
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One motivation for the analysis in [Kl, E^ was to test the large- A^c prediction 



Lg = 2Lp The above result shows that a small non-zero value is preferred. To obtain 
a clean error analysis, we have repeated the fitting procedure using the variables 

Xi = L2 — 2L^ — L3 , 

Xr T r 
2 ~ -^2 ? 

X3 = {Ll-2L[)/L3 . (5.82) 

We performed a fit to Ke4 and tttt data, including the theoretical error in G as 
discussed above, and found 

Xi = (4.8 ±0.8) -10"^ , 

X3 = -0.17+°i2 . (5.83) 

The result is that the large- Xc prediction works remarkably well. 



5.9 Predictions 

In this section we make several predictions using the L['s from table 1 in Ref. 0. It 
is clear that new and more accurate data on ii"e4 will allow for a better determination 
of L\, L2 and L3, and may correspondingly modify our predictions. However, unless 
a dramatic change in the values of these constants occurs, the modified predictions 
will be within the errors that we give. 

Whereas the slope A^ was assumed to coincide with the slope Xf in the final 
analysis of the data in Ref. [Q , these two quantities may differ in the chiral repre- 
sentation. Furthermore, our amplitudes allow us to evaluate partial and total decay 
rates. In this section, we consider the slope A^ and the total rates. 



The slope A^ 



We consider the form factor g introduced in (|5.78| ) and determine its slope A 



'9 



-g{s^, si) = g{4M^, sl){l + \,{si)q' + 0(g^)) (5.84) 

from the one-loop expression for G. The result is Ag(0) = 0.08. As the slope is a 
one- loop effect, higher-order corrections may affect its value substantially. For this 
reason, we have evaluated \g also from the modified form factor obtained by using 
the complete resonance propagators (and the corresponding Lj's), see Ref. |^^. The 
change is AA^ = 0.025. We believe this to be a generous error estimate and obtain 
in this manner 

Ag(0) = 0.08 ±0.025 . (5.85) 

The central value indeed agrees with the slope A in (^.37|) . 
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Table 5.3: Approximations used to evaluate the total rates in table |5.4| . Use of 
/ = /cHPT, 9 = ^CHPT reproduces the one-loop results in table E^ to about 1%. 



a) K^ decays 





vrOvrOe+z/e 


TT+TT /U+Z/^ 


ttOttV+z^m 


F2 

Fs 
F, 


-Xf 




{PL)f 


Xf + a^{PL) COS e^g 



- { (PL) / + s;/?cHPT + cr^Xcose^g] 


-Xf 






b) K^ decays. Shown are the amplitudes divided by v2- 





7r°7r e+z/e 


7r°7r 11+ u^ 


Fi 
F2 
F3 

F, 


^^c~HPT + (^TriPL) cosO^g 

(T^{s^siY'^g 



-{{PL)F-^^^ + a^Xcose^-g} 


XF-jjpT + a^iPL) cos 9^g 

(Tnis^siY/^g 



- {(PL)F-^p^ + siR-^^^ + or^Xcos^,^} 



Total rates 

Once the leading partial waves / and g are known from e.g. K^ -^ -K^-K^e^Ve decays, 
the chiral representation allows one to predict the remaining rates within rather 
small uncertainties. We illustrate the procedure for K^ -^ Ti^Ti^e^Ue. According 
to Eq. ( |5.24| ), the relevant amplitude is determined by F^,G^,R'^ and H~. The 
contribution from H is kinematically strongly suppressed and completely negligible 
in all total rates, whereas the contribution from R is negligible in the electron modes. 
Using the chiral representation of the amplitudes F~^ and G~, we find that the rate 
is reproduced to about 1%, if one neglects G~ altogether and uses only the leading 
partial wave in the remaining amplitude, F^ ~ —Xf. From the measured [|53| form 
factor / = 5.59(1 + 0.08g^) we then find TK+^TTO-n-Oe+iye = 1625sec~^. Finally, we 
estimate the error from 



AF = {[F(/,(0) + A/„A;)-F(/,(0),A;)f+ 



1/2 



[F(/,(0), \f + A\f) - F(/,(0), \f)]'} = 90sec-i 



(5.86) 



where A/s = 0.14, AA/ = 0.02. The final result for the rate is shown in the row 
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Table 5.4: Total decay rates in sec ^. To evaluate the rates at one-loop accuracy, we 
have used the L['s from table 1 in Ref. |0. The final predictions are evaluated with 
the amplitudes shown in table U' using / = 5.59(1 ± 0.08q'^),g = 4.77(1 ± O.OSq^). 
For the evaluation of the uncertainties in the rates see text. 



a) K^ decays 





TT + TT e+Z/g 


n^TT^e+Ue 


TT + TT^/i+Z/^ 


TT^'fl+U, 


tree 
one-loop 


1297 
2447 


683 
1301 


155 

288 


102 
189 


final 
prediction 


input 


1625 
±90 


333 

±15 


225 
±11 


experiment 

i 


3160 
±140 


1700 
±320 


1130 
±730 





b) K^ decays 





ir^n e+z/g tt'^tt fJ^'^i'fj. 


tree 
one-loop 


561 55 
953 94 


final 
prediction 


917 88 
±170 ±22 


experiment 


998 
±39 ± 43 



"final prediction" in table |5.4| , where we have also listed the tree and the one-loop 
result, together with the experimental data. The evaluation of the remaining rates 
is done in a similar manner - see table ^]3| for the simplifications used and table |5]4 
for the corresponding predictions. 

We have assessed an uncertainty due to contributions from F~^prj,, -Rchpt in the 
following manner, i) We have checked that the results barely change by using the 
tree level expression for -Rchpt instead of its one-loop representation. We conclude 
from this that the uncertainties in -Rchpt do not matter, ii) The uncertainty from 
F~^p^ is taken into account by adding to AF in quadrature the change obtained by 
evaluating -F^pt ^^ -^3 = —3.5 ± 1.1 = —2.4. iii) In K^ decays, we have also added 



in quadrature the difference generated by evaluating the rate with M^r 

The decay K^ — >■ n^TT'e^Ue has recently been measured | 
higher statistics than before P] . We display the result for the rate in the first column 



135 MeV. 
with considerably 



of table |5.4| b. The quoted errors correspond to the errors in the branching ratio ||6^ 
and do not include the uncertainty in the total decay rate quoted by the PDG B. 
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Notice that the value for L3 determined in |pO[ should be multiplied with — 1 Wt 



5.10 Improvements at DA$NE 

The chiral analysis of Ki^ decays has been used so far for three purposes: 



1. The KeA data from Ref. [Q allows one to make predictions for the slope 
of the G form factor, for the total rates in all the channels and for the vrvr 
scattering lengths. These are given in Eq. ( 5.85| ), in table ^.41 and in table 5.2 , 
respectively. 



2. The same K^a data allow one to test the large- iVc prediction, see Eqs. ( |5.82| ) 
and ( ^1831) ■ 



3. The full set of Kei and tttt scattering data allows the best determination of 
the coefficients L\,L2 and L3 in the chiral Lagrangian, see ( ^.80|) . 



In the next generation of Ki^ decay experiments, there is the opportunity to 
improve the phenomenology of Ki^ (see table |5.1|): 



1. A very useful innovation would be to analyze the experimental data with a 
modified chiral representation. In the latter, the full S- and P- wave parts of Fi 
and F2 could be inserted, using the chiral representation solely to describe the 
small background effects due to higher partial waves / > 2. To be more precise, 
one would take for R and H the one-loop chiral representation, whereas for G 
one writes 

G = g{s^,si)e''- + ^G , 
AG = GcHPT-— /c/f^sin^^^GcHPT , (5.87) 

and similarly for F. The unknown amplitudes g(sT,,si), fs{sT,,si) and the 
phases 6p, 6s would then be determined from the data. We have checked that, 
if the errors in the form factors determined in this manner can be reduced by 
e.g. a factor 3 with respect to the ones shown in ( |5.37| ), one could pin down 
particular combinations of L^, Lg and L3 to considerably better precision than 
was shown above. This is true independently of an eventual improvement in 
the theoretical determination of the higher-order corrections in the form factor 
G - which is a theoretical challenge in any case. 

2. The present experimental uncertainty on G is still too large to provide a precise 
value for the large- iVc parameter (Lg — 2L\)/L^. {K^ — >■ ir^TT'e^Ue decays are 
mainly sensitive to G^_ which in turn can be used to pin down L3. K^ — >■ 
7r°7r°e"'"z/e is mainly sensitive to F^_ which contains Li, L2 and L3.) 
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3. The observation of all Ku reactions with high statistics could provide a cleaner 
separation of the various isospin amplitudes. 

4. Finally, we come to a most important point. As we mentioned already, 



K^ -^ TT+TT e+z/g has been used [59| to determine the isoscalar S-wave scat- 



tering length with the result a^ = 0.26 ± 0.05. This value must be compared 
with the SU(2)xSU(2) prediction (T^, |78[ a[] = 0.20 ± 0.01. Low-energy vrvr 



scattering is one of the few places where chiral symmetry allows one to make 
a precise prediction within the framework of QCD. In their article, Rosselet et 
al. 

79| Oq = 0.16 in the following manner: 



comment about the discrepancy between Oq = 0.26 ±0.05 and the leading- 



order result 



it appears that this 
prediction can be revised without any fundamental change in current algebra 
or in the partial conservation of axial- vector current [BO, Bl[." Today, we know 



that this is not the case: The standard picture of the vacuum structure in 
QCD ||S^ would have to be revised, should the central value Oq = 0.26 be 
confirmed with a substantially smaller error. For recent work which supports 
this scenario see the contribution of Knecht and Stern in this Handbook 



Ki4 decays are - at present W^ - the only available source of clean information 



on TTTT S'-wave scattering near threshold. We refer the reader to Ref. |8J] for a 
detailed analysis of the issue. 
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6 i^eS decays 

In this subsection we discuss the decays 



K+ - 


-* TT^TT n^e'^Ue 


K+ - 


-^ 7r\%'>e+u, 


K' - 


-» 7r°7r°7r"e^z/e 


K' - 


-» n^n^'K^e'^h'f 



(6.1) 

We do not consider isospin violating contributions and correspondingly set niu = md, 
a = 0. 



6.1 Matrix elements and decay rates 

The matrix element for K -^ TiTme^Ue is 

T = ^v:MPuh,ii - iMPeW^' - An , (6.2) 



where 



/^ =< nip,)7rip,)nip,)ont\C-''iO)\K{p) > ; I = V, A. (6.3) 



The decay rate is calculated from 

d'^ = 2M f27r)" ^ \T\^dLIPs{P]Pe,Pu,Pl,P2,P3) ■ 
^ V / spins 

6.2 Previous experiments 

The Particle Data Group ^ quotes the upper bound 

BR{K+ -^ 7r°7r°7r°e+z/e) < 3.5 • 10"^ 

6.3 Theory 

In CHPT, the leading-order contribution is given by the matrix element of the vector 



current. The corresponding rates are displayed in table (|6.1| ). The smallness of these 
rates is due to the suppression of phase space. Indeed, consider the ratio of the four- 
and five- dimensional phase space volumes in the neutral pion channel, 

MJ,JdLisp{p;Pe,P.,Pi,P2)/2\ (27r)i^ ~ 2 3 • 10^ 

(27r)12 I dLISpiP^,Pe,Pu,Pl,P2,P3)/^^■ 

It agrees well with the ratio of the corresponding rates at tree level. 



^The material in this section is taken from Ref. [p5[ . 
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Table 6.1: Rates of i^es decays, evaluated from the leading-order term. 





branching ratio 


K+ - 


-* TT + TT 7r°e+Z/e 


3-10-12 


K+- 


-^ TrOTrOvrOe+z/e 


2.5 ■ 10-12 


K'- 


->■ 71^71^71^ e^Ue 


12 ■ 10-12 


K'- 


->■ TT + TT^TT^e + Z/e 


33 ■ 10-12 



^(ir^ 



7r°7r°e+Z/e)iree 



^(K^ 



^Oy^Oy^Og+w 



3.4-10' 



e I tree 



(The corresponding ratios for K^ -^ n^e'^Ue/K^ -^ n^Tc'^e^Ve are 1.4 ■ 10^ and 
0.53 ■ 10^ for phase space volumes and decay rates, respectively.) The contributions 
at order p^ are due to i) the corrections to the matrix element of the vector cur- 
rent, and to ii) the matrix element of the axial current. The latter stems from the 
Wess-Zumino-Witten Lagrangian Cwzw- Both the local and nonlocal term in the 
anomalous action contribute. The nonlocal part is suppressed by the factor irie in 
the matrix element (in addition to the phase space suppression just mentioned). 

Based on our experience with KeA decays, we expect the terms of order p^ to 
enhance the above rates by roughly a factor of two to three. 



6.4 Improvements at DA$NE 

According to the standard model, K^.^ decays are invisible at DA$NE, but the 
existing upper limits can be improved significantly. 
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A Notation 

The notation for phase space is the one without the factors of 27r. For the decay 
rate of a particle with four momentum p into n particles with momenta pi, . . . ,Pn 
this is 

dLiPsip;pi, ...,Pn) = 5^ip - II Pi) n -rri- ■ (^-1) 

We use a covariant normalization of one-particle states, 

<p'\p>={27r)^2p'^6^{p' -p) , (A.2) 

together with the spinor normalization 

u{p,r)u{p,s) = 2m6rs ■ (A. 3) 

The kinematical function X{x, y, z) is defined as 

A(x, y, z) = x^ + y^ + z^ — 2{xy + yz + zx) . (A. 4) 

We take the standard model in the current x current form, i.e., we neglect the 
momentum dependence of the IV-propagator. The currents used in the text are : 



T/4-J5 


= <lli^id^i-i^^)<l = sii^u 


/l4-i5 


= 97^752(^4 -«A5)g = S7^75« 


-irem 


= ql^lQ(l 


Q 


= diag(2/3,-l/3,-l/3). 



(A.5) 

The numerical values used in the programs are the physical masses for the particles 
as given by the Particle Data Group 0. In addition we have used the values for 
the decay constants derived from the most recent measured charged pion and kaon 
semileptonic decay rates [|], ^ : 

F^ = 93.2 MeV 

Fk = 113.6 MeV. (A.6) 

We do not need values for the quark masses. For the processes considered in this 
report we can always use the lowest order relations to rewrite them in terms of the 
pseudoscalar meson masses (see Ref. 0). For the KM matrix element we used the 
central value \Vus\ = 0.220 of Ref. 0. The numerical values for the L^-{Mp) are 
those given in table 1 in Ref. 0. 

The number of events quoted for DA$NE are based on a luminosity of 5 ■ 
10'^^ cm^^s^^, which is equivalent |]1|] to an annual rate of 9 ■ 10^ (1.1 ■ 10^) tagged 
K"^ {Kl) (1 year = 10'^ s assumed). 
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Whenever we quote a branching ratio for a semileptonic K^ decay, it stands for 
the branching ratio of the corresponding K^ decay, e.g., 

BR{K^ -^ TT-l+u) = BR{Kl -^ n^l^u) . (A.7) 

We use the Condon-Shortley phase conventions throughout. 
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B Loop integrals 

In this appendix we define the functions appearing in the loop integrals used in 
the text. First we define the functions needed for loops with two propagators, 



mainly in the form given in Ref. |^. We consider a loop with two masses, M 
and m. All needed functions can be given in terms of the subtracted scalar integral 
J{t) = Jit) - J(0), 

/dp 1 

with t = k'^. The functions used in the text are then : 

J{t) = o / dx log — ^ -^ 

^ ^ 167r2 Jo ^ M2 - Ax 

1 [^ A, m2 E, m^ v^ , (t + VX)^-A^] 
^ 2 H log — -r — — log — log ' 



327r2 [ t ^ M2 A '^ M2 t "^ [t- v^)2 - A2 

r{t) = J{t) - 2k , 



12t ' ' ' ' 3t2 ' ' 2887r2 6 

^ M'^rn? , m2 1 
S + 2^^ log ■ 



967r2t [ A "^ M2 J ' 

A2 - 

ir(t) = lj(t), 

H{t) = \^,t + ^^[tM^{t)-L{t)l 

A = M"^ ~m^ , 
S = M^ + m^ , 
A = A(t,M2,m2) = (t + A)2-4tM2. (B.2) 

In the text these are used with subscripts, 

J^.[t) = J{t) with M = Mi,m = Mj (B.3) 

and similarly for the other symbols. The subtraction point dependent part is con- 
tained in the constant k 

1 M2log(^)-m2log(4) 
where /i is the subtraction scale. 
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In addition, in subsection ^ these functions and symbols appear in a summation 
over loops / with 

Jj(t) = J{t) with M = Mj,m = mj ; 
E/ = M'j + m] (B.5) 

and again similarly for the others. There the combination B2 appears as well : 

1 , , A(t,M2,m2)J(t) tS-SM^m^ M^ 

(3S-t) j^^ + — ^^^_^^ — log^ry • 



2887r2 ' ' I2t SSAtt^A ^ m^ 

The last formula to be defined is the three propagator loop integral function 
C{ti, t2, M'^,m'^) where one of the three external momenta has zero mass and two of 
the propagators have the same mass M. Here ti = (gi + ^2)^, t2 = g| and qf = 0. 

2 2n , f d^P 1 



C{ti,t2,M\m 



{27tY (p2 - M2)((p + gi)2 - M2)((p + gi + 52)^ - m^ 

^1 rl—x 

dx dy- 

Jo Jo 



167r2 Jo Jo ^ M2 - y{A + ti) + xy{ti - ts) + 2/^ti 
-Lz2 — T—r - Li2 



\y+iti)J Ky-ih)/} 

y±{t) = ^{t + A±^A(t,M2,m2)| (B.6) 

where Li2 is the dilogarithm 

Lt2ix) = - f'^\ogil-xy) . (B.7) 

Jo y 
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C Decomposition of the hadronic tensors I^^ 

Here we consider the tensors 

J/^- = J dxe^'^^+^'^y < I TVUx)IL^M I K+ip) > , I = V,A (C.l) 



and detail their connection with the matrix element ( |1.2|) . 

The general decomposition of A^^", V^" in terms of Lorentz invariant amplitudes 



reads P, 11| for g^ 7^ 



1=V^'' = iVie^^'^^q^pp (C.3) 

where the form factors Ai^q^^W"^) and Vi{q^,W'^) are analytic functions of q^ and 
W^. Fy{q^) denotes the electromagnetic form factor of the kaon (F^(0) = 1). A^^ 
satisfies the Ward identity 

q^A^- = -V2FKP''. (C.4) 



and 



In the process ( |1 . 1| ) the photon is real. As a consequence of this, only the two 
form factors Ai{0, W'^) and Vi(0, W'^) contribute. We set 

A{W^) = Ai{Q,W^) 

V{W^) = Vi{0,W^) (C.5) 



and obtain for the matrix element (|1.2| ) 

T = -zGp/V2eVus*e; {V2FkI'1 - (y^^ - A^L}^ , (C.6) 

with 

n = /" + m,tx(p.)(l+75) ?'^^]_\A pi) ■ (C-7) 

^i - \Pi + ?) 

Grouping terms into an IB and a SD piece gives ( |1.2| , |1.3| ). As a consequence of 
( |C.4| ), T is invariant under the gauge transformation e^ — > e^ + g^. 

The amplitudes Ai, A2 and Vi are related to the corresponding quantities Fa, R 
and Fv used by the PDG g by 

- V2Mk{Ai, A2, V,) = {Fa, R, Fy). (C.8) 
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The last term in ( |C2| ) is omitted in p. It contributes to processes with a virtual 
photon, K^ -^ l^vil'+l'-. 

Finally, the relation to the notation used in [^ Q is 

2{A±Vf = {a,±Vkf [4] 
V2{A,V) = (Fa^Fv) § . (C.9) 
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D One-loop corrections to K14 form factors 

In this appendix we give the expression of the unitarity corrections to the form 
factors F, G, and i? (13,11, |K 



UF{s^,t,u) 



Ao{s^) + AF{t) + B{t,u) 



(D.l) 



with 
Apit) 



1 Ss M^ 



+ 



16 
1 



4 



(14Mi + UMl - m)J],^{t) + (2Mi + 2Ml - 3t)j:^At) 



{ml - IMl + ^t)KK.{t) + {Ml - ^Ml + 2,t)Kr,K{t) 



9{LKAt) + L^K{t)) + (3Mi - ?,Ml - 9t)(M^^(t) + M„V(t)) 



B{t,u) 



-l^{Ml + Ml-t)rj,,{t)-{t 



u 



(D.2) 



The loop integrals J^^(s7r), . . . which occur in these expressions are listed in appendix 
1^. The functions Jpq and Mpg depend on the scale fi at which the loops are 
renormalized. The scale drops out in the expression for the full amplitude. 

The imaginary part of F^^Ao(s7r) contains the / = 0, S-wave tttt phase shift 



6^o{s^) = {327rF^)-'{2s^ - M^)a^ + 0{E' 



(D.3) 



as well as contributions from KK and rjrj intermediate states. The functions AF{t) 
and B{t, u) are real in the physical region. 



UG{s^,t,u) 



Ai{s^) + AG{t) + B{t,u) 



(D.4) 



with 



AifsJ 



2s^ { M:,{s^) + -MW{s.) 



AG{t) = -[{2Ml + 2Ml + ?>t).r^^{t)-{2Ml + 2Ml-3t)Jl,^{t) 



+ 



:-3Mi + 7Ml - 5t)KKn{t) + {-Ml + SM^ - 3t)K,K{t) 



LK.{t) + L,K{t) - {Ml - Ml + t){M],^{t) + M:^j,{t)) 



(D.5) 



The imaginary part of F^ ^Ai(s7r) contains the / = 1, P-wave phase shift 

5\{s^) = {967rFl)-'s^al + 0{E^) . (D.6) 
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as well as contributions from KK intermediate states. The function Aq is real in 
the physical region. 

The unitarity corrections Uz-, Uq in the form factor R in ( [5.58| ) are 



U 



z = s^Ao(s^)+z/Ai(s^)--M|m2j; (s^) 

y 



^ 32 



ll(s^-z/)2-20S(s^-z/) 



12S^ 



■JK.it) 



+ l[3(.^-z.)-2S]^j;^(t) 



1 

4' 



+ -{s^ + iyyr^^{u) 






+ 



2s^iiy + 4S) - 3si + u' ~ IQM^M'^ M],^(t) + M;^(t) 



- -{'is^ + v-2'L){L^K{t) + LK.it)) 



(D.7) 



- l{Hsn-i^) + HMl-si)-6J:)KKAt) 



- -(3(s,-z/) + 3(M^-sO-2S)ir,i,(t) 

o 

9 

- -{L^K{t) + LxAt)) 

+ l{A{u + 2Ml) - ?,{Ml - s,)){MUt) + M;^(t)) 



(D.8) 



with 



E = M^ + M^ . 
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